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SECTION  1 
INTRODUCTION 


Many  different  military  targets  must  be  analyzed  with  regard  to 
their  vulnerability  to  the  blast  effects  from  a  nuclear  weapon.  Often 
these  targets  are  comprised  of  structural  elements  such  as  plates, 
beams,  and  cylindrical  bodies.  A  convenient,  simple,  and  approximate 
means  of  defining  the  vulnerability  of  a  structural  element  is  by 
Pressure  (P)  -  Impulse  (I)  curves  for  different  levels  of  damage.  In 
the  past  these  P-I  curves  have  been  generated  experimentally  at  con¬ 
siderable  cost.  The  same  sort  of  thing  can  be  accomplished  by  using 
complex  computer  codes.  Consequently,  the  basic  objectives  of  the 
current  study  are  the  following:  1)  to  develop  P-t  curves  from  complex 
computer  codes  for  different  structural  shapes  such  as  beams,  plates, 
and  cylinders  for  three  levels  of  damage  -  light,  moderate,  and  severe, 
and  2)  to  include  these  P-I  curves  in  a  computer  code  which  the  analyst 
can  use  to  rapidly  approximate  the  vulnerability  of  a  specific  target  to 
blast. 

The  response  of  structural  systems  to  nuclear  blast  effects  has 
been  of  interest  to  the  military  community  for  many  years  and  past 
efforts  aimed  at  assessing  structural  blast  damage  have  involved  both 
experimental  and  analytical  approaches.  Since  full-scale  structural 
blast  damage  response  cannot  be  measured  until  the  structural  system  has 
been  built  and  experiments  performed,  experimental  studies  tend  to  be 
expensive,  particularly  if  tests  of  a  variety  of  structural  configura¬ 
tions  and  damage  levels  are  required.  On  the  other  hand,  analytical 
assessment  of  structural  blast  damage  can  be  accomplished  either  after 
the  actual  structure  has  been  built  or  before  production,  while  the 
structure  is  still  in  the  design  stage.  Thus,  it  is  evident  that  the 
availability  of  analytical  data  for  different  damage  levels  of  struc¬ 
tural  blast  response  on  a  variety  of  structural  configurations  would 
certainly  be  of  interest  to  both  weapon  systems  designers  and  analysts, 
particularly  for  quick-look  vulnerability  assessment.  Consequently,  the 
U.S.  Army  Ballistic  Research  Laboratories  sponsored  the  present  analy¬ 
tical  blast  study  program  whose  ultimate  objective  was  to  make  such  data 
readily  available. 

The  basic  analytical  technique  which  was  selected  to  provide  the 
required  blast  damage  assessment  data  was  the  pressure-impulse  approach 
which  had  proven  quite  successful  in  previous  structural  blast  damage 
assessment.  Section  2  of  this  report  provides  a  background  discussion 
of  the  pressure- impulse  approach  and  its  application  to  the  proposed 
program  objective.  As  specified  in  the  contract  statement  of  work,  the 
first  detailed  requirement  of  the  program  was  to  review  selected  mili¬ 
tary  structural  configurations  and  to  select  typical  simple  structural 
shapes  which  could  be  interpreted  as  representative  component  shapes 

11 

^  phbckdino  PAcat  blank-not  filmed 


for  those  systems.  Once  the  representative  component  shapes  were  cho¬ 
sen,  the  next  step  specified  was  to  determine  what  response/damage 
levels,  materials,  boundary  conditions,  and  loading  functions  would  be 
of  interest  in  evaluating  blast  effects  on  those  representative  compo¬ 
nents  and  to  utilize  appropriate  computer  codes  in  conjunction  with  the 
pressure-impulse  approach  to  calculate  the  required  response  data.  In 
Section  3  of  the  report,  descriptions  are  given  of  the  type,  scope  and 
limitations  of  the  pressure-impulse  analyses  which  were  undertaken. 

These  descriptions  include,  among  other  things,  the  types  of  and  con¬ 
straints  on  the  structural  component  models  considered,  the  applica¬ 
bility  of  the  assumed  loading  functions,  the  computer  codes  employed  in 
the  structural  response  investigations,  and  the  implementation  of  the 
basic  pressure-impulse  approach. 

Section  4  presents  an  explanation  of  the  manner  in  which  analytical 
fitting  functions  were  used  to  process  the  derived  pressure-impulse  data 
into  a  form  suitable  for  incorporation  into  the  vulnerability  computer 
code  which  was  formulated  as  the  specified  end  requirement  of  the 
contract.  A  description  of  the  code  together  with  an  explanation  of  its 
usage  for  rapid  assessment  of  blast  structural  vulnerability  is  pre¬ 
sented  in  Section  5.  Finally,  Section  6  presents  the  conclusions  and 
recommendations . 
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SECTION  2 

PRESSURE- IMPULSE  APPROACH  TO  STRUCTURAL  BUST 
DAMAGE  DEFINITION 


The  response  of  military  structures  to  nuclear  weapon  effects  has 
been  of  interest  ever  since  the  first  atomic  weapon  was  introduced  into 
the  offensive  weapon  arsenal.  In  particular,  military  strategists  have 
been  faced  with  the  problem  of  assessing  the  size  of  nuclear  weapons 
needed  to  destroy  the  enemy's  offensive  and  defensive  capabilities  and, 
conversely,  how  to  assess  the  hardness  of  one's  own  military  structures 
to  nuclear  weapons  effects. 

The  work  reported  herein  describes  one  study  effort  performed  to 
improve  the  capability  for  rapid  prediction  of  damage  to  mili¬ 
tary  equipment  caused  by  the  blast  effects  resulting  from  nuclear 
weapons  detonations.  More  specifically,  the  intent  of  the  study  is  to 
provide  nuclear  blast  effects  methodology  which  can  be  utilized  in 
either  preliminary  design  of  structural  systems  for  blast  effects  or 
in  making  a  rapid  assessment  of  nuclear  blast  damage  to  existing  struc¬ 
tural  systems. 

The  basic  concept  of  blast  effects  damage  assessment  employed  in 
the  present  report  is  the  pressure-impulse  approach  for  relating  dynamic 
structural  response /damage  to  structural  blast  loadings  created  by 
nuclear  weapons  detonations.  As  will  be  subsequently  discussed,  the 
pressure-impulse  approach,  hereafter  referred  to  as  the  P-I  approach, 
is  not  a  new  concept  in  structural  damage  assessment  but  has  been  used 
for  many  years  by  structural  analysts  interested  in  structural  blast 
damage  assessment. 

In  structural  analysis,  the  response  of  a  particular  structural 
item  is  functionally  dependent  on  many  parameters,  including  the  geome¬ 
try  of  the  structural  item,  the  material  used,  the  structural  support 
fixity,  and  the  nature  of  the  loading.  If  the  load  imposed  on  the 
structure  does  not  vary  with  time,  the  structure  may  be  analyzed  as  a 
static  structure  and  handbook  formulas  are  often  available  for  deter¬ 
mining  the  structural  deformation  as  a  function  of  the  magnitude  and 
spatial  distribution  of  the  applied  load.  On  the  other  hand,  if  the 
load  imposed  on  the  structure  is  time-varying  then  dynamic  structural 
analysis  is  required  since  the  response  is  functionally  dependent  on 
the  temporal  variation  of  the  loading  as  well  as  the  dynamic  charac¬ 
teristics  of  the  structure.  It  is  the  latter  case  which  is  of  interest 
in  the  present  report  since  the  blast-type  loadings  which  are  of 
interest  in  this  program  are  dynamic  loadings  rather  than  static  load¬ 
ings. 
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The  detonation  of  an  explosive  charge  creates  a  shock  front  or 
blast  wave  which  propogates  radially  outward  from  the  center  of  the 
detonation.  It  is  not  the  intent  of  this  report  to  present  a  detailed 
explanation  of  shock-wave  propagation  and  it  suffices  to  say  that  at  the 
advancing  shock  front  there  is  an  almost  instantaneous  increase  in  air 
pressure  above  the  ambient  air  pressure,  called  the  overpressure. 

As  the  distance  from  the  detonation  center  increases,  the  peak 
overpressure  decreases.  At  a  given  distance,  the  magnitude  of  the 
overpressure  behind  the  shock  front  decreases  with  increasing  time  after 
detonation,  so  that  the  time-history  of  the  pressure  experienced  at  a 
point  might  appear  as  depicted  in  Figure  2.1,  where  the  overpressure  is 
plotted  versus  the  time  after  detonation.  The  time  t^  represents  the 
time  required  for  the  shock  front  to  travel  from  the  detonation  center 
to  the  point  in  question,  the  time  t2  represents  the  positive  phase  of 
the  pressure  pulse  (i.e.,  the  time  duration  for  which  the  pressure  is 
greater  than  ambient  pressure),  and  time  tj  represents  the  negative 
phase  of  the  pressure  pulse  during  which  time  the  pressure  is  less  than 
ambient  pressure.  In  this  report,  only  the  positive  phase  of  the 
overpressure-time  history  is  of  significance  and,  hence,  only  blast 
loadings  corresponding  to  the  time  period  t2  in  Figure  2.1  are  considered. 
The  blast  loadings  of  interest,  then,  are  seen  to  be  pressure  loadings 
which  change  over  a  finite  time  period  and  are  dynamic  rather  than 
static  in  character. 


If  the  typical  positive  phase  pressure  loading  depicted  in  Figure  2.1 
is  applied  to  a  particular  structural  item,  such  as  a  fixed-end  beam  for 
instance,  the  deformation  pattern  of  the  beam  will  vary  with  time  since, 
as  previously  indicated,  the  dynamic  response  is  a  function  of  both  the 
dynamic  characteristics  of  the  structure  and  the  time-history  of  the 
loading  function.  Within  some  finite  period  of  the  structural  response, 
however,  maximum  values  of  deflections,  strains  and  stresses  will  be 
reached  and  the  maximum  amount  of  damage  which  is  incurred  on  the  struc¬ 
ture  can  be  related  to  one  of  the  maximum  response  quantities.  If  the 
same  structure  is  subjected  to  another  pressure  pulse  of  the  form 
depicted  in  Figure  2.1  but  with  a  di  .erent  magnitude  of  either  peak 
overpressure  or  positive  duration,  the  maximum  value  of  the  same  response 
parameter  will  differ  from  the  peak  value  experienced  during  the  original 
loading.  There  are,  however,  other  combinations  of  peak  overpressure 
and  positive  duration  that  will  result  in  the  same  peak  response  or  peak 
damage  level  as  that  obtained  from  the  original  loading. 

Referring  back  to  the  loading  depicted  in  Figure  2.1,  to  be  more 

consistent  with  the  literature,  Ap  will  henceforth  be  referred  to  as 

max 

Apg,  the  value  of  overpressure  at  the  shock  front,  and  t 2  will  be 

referred  to  as  t0,  the  positive  phase  duration.  Since,  as  shown  in 
Figure  2.1,  the  overpressure  at  a  point  removed  from  the  detonation 
center  decays  in  an  exponential  manner  with  time,  an  expression  for  the 
overpressure  as  a  function  of  time  may  be  written  as 


Ap(t) 


(2.1) 
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FIGURE  2.1 


where  t=0  corresponds  to  the  shock  front  arrival  at  the  point.  It  can 
be  seen  then  that  the  entire  loading  time-history  curve  can  be  defined 
by  only  two  parameters,  Aps  and  tQ.  It  would  be  entirely  possible  then 
to  load  the  identical  structure  with  the  same  forcing  function  of  Eq. 
(2.1),  but  with  varying  combinations  of  Aps  and  tQ,  and  determine  the 
combinations  of  the  two  parameters  which  would  result  in  the  same 
response  level.  The  plot  of  these  various  Aps  and  tD  pairs  would  then 
define  a  curve  of  equivalent  level  of  response  or  damage,  commonly 
referred  to  as  an  iso-damage  curve. 

For  forcing  functions  of  a  more  general  shape  than  that  expressed 
by  Eq.  (2.1),  a  greater  number  of  parameters  are  required  to  define  the 
loading  time-history.  Noting  that  the  response  will  depend  upon  any 
parameter  which  affects  the  loading  time-history  (up  to  the  time  of 
maximum  response),  it  becomes  apparent  that  the  response  will  then 
depend  upon  more  than  two  parameters  for  an  arbitrary  loading  time- 
history.  The  construction  of  iso-damage  curves  depending  upon  more  than 
two  parameters  becomes  cumbersome,  however.  Efforts  in  the  P-I  area 
have  been  directed  toward  selecting  the  two  most  dominant  loading 
parameters  in  the  hope  that  the  influence  of  the  unused  parameters  will 
be  small. 

The  most  fundamental  approach  originated  in  the  work  of  Sperrazza 
(Ref.  1)  in  1951.  Sperrazza  based  his  iso-damage  curves  on  the  two  key 
parameters  of  peak  pressure  of  the  forcing  function  and  the  total 
impulse  associated  with  the  forcing  function  time  history,  hence  the 
resulting  name  pressure-impulse  (P-I)  approach.  In  a  later  study  repor¬ 
ted  in  Ref.  2,  Youngdahl  noted  that  the  structural  response  was  strongly 
dependent  on  the  shape  of  the  forcing  function  even  for  pulses  which 
have  the  same  total  impulse  and  peak  pressure,  which  is  implied  by  the 
foregoing  comments.  In  his  paper,  Youngdahl  suggested  the  use  of  an 
"effective  pressure"  and  an  "effective  impulse",  showing  that,  in  some 
sense,  these  two  parameters  represented  a  better  choice  than  peak 
pressure  and  fatal  impulse.  Certain  procedural  difficulties  occur  in 
employing  Youngdahl' s  definitions  of  effective  pressure  and  effective 
impulse.  Hence,  in  a  later  study.  Ref.  3,  Cummings  and  Schumacher 
suggest  a  simplified  version  of  Youngdahl 's  approach.  In  which  the 
"effective  pressure"  and  "effective  impulse"  are  redefined  by  changing 
the  integration  limits  involved.  This  approach,  along  with  the  funda¬ 
mental  P-I  approach  employing  peak  pressure  and  total  impulse,  is  used 
in  the  present  study. 


Ref.  1.  Sperrazza,  J.,  Dependence  of  External  Blast  Damage  to  the  A-25 
Aircraft  on  Peak  Pressure  and  Impulse,  Memorandum  Report  575, 

USA  Ballistic  Research  Laboratories,  APG,  Maryland,  Segtember^l951. 

Ref.  2.  Youngdahl,  C.K.,  "Correlation  Parameters  for  Eliminating  the 

Effect  of  Pulse  Shape  on  Dynamic  Plastic  Deformation",  Journal  of 
Applied  Mechanics,  Transactions  of  the  ASME,  pp.  744-752, 

September  1970. 

Ref.  3.  Schumacher,  R.,  and  Cummings,  B.,  A  Modified  Pressure-Impulse  Blast 
Damage  Model,  BRL  MR  2724,  AD  No.  A036196,  USA  Ballistic  Research 
Laboratory,  APG,  Maryland,  January  1977. 


16 


It  is  important  to  recognize  that  any  response  which  is  found  as  a 
function  of  only  two  parameters  of  an  arbitrary  loading  can  be  only  an 
approximation  to  the  true  response.  If  iso-damage  curves  calculated 
based  upon  the  loading  time-history  of  Eq.  (2.1)  are  applied  for  loading 
time-histories  with  similar  characteristics,  the  approximation  will  be 
good.  If  they  were  applied,  for  example,  to  a  half  sine-wave  forcing 
function,  the  approximation  would  be  expected  to  be  poor.  The  data 
generated  are  based  upon  the  loading  time-history  of  Eq.  (2.1),  and  must 
therefore  be  applied  only  to  similar  loadings.  Since  blast  loadings, 
even  including  reflection-diffraction  effects,  are  basically  similar  to 
that  given  by  Eq.  (2.1),  the  approximation  involved  in  the  present 
context  in  ignoring  the  complete  details  of  the  loading  time-history 
should  be  acceptable.  Appendix  A  shows  the  type  of  variation  that  may 
be  expected  in  the  P-I  curve  for  different  forcing  functions  within 
the  same  class. 
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SECTION  3 

DEVELOPMENT  OF  P-I  CURVES 


3.1  Introduction 

As  indicated  in  Section  1,  the  primary  objectives  of  the  technical 
effort  reported  herein  were  to  develop  P-I  curves  for  simple  structural 
elements  subjected  to  nuclear  blast  overpressure  loadings  and  to  include 
these  P-I  curves  in  a  vulnerability  assessment  code.  A  description  of 
the  structural  elements  selected  for  analysis,  the  modeling  of  the  ele¬ 
ments  including  the  boundary  conditions  and  materials  assumed,  the 
assumed  aerodynamic  loading  functions,  the  selection  of  the  response 
damage  criteria,  and  the  methods  of  response  analysis  are  discussed  in 
the  first  portion  of  this  section.  The  remainder  of  the  section  is 
primarily  devoted  to  the  formulation  of  the  P-I  curves  which  were  needed 
for  inclusion  in  the  vulnerability  code. 

3.2  Selection  of  Typical  Structural  Models 

The  process  of  selecting  representative  structural  items  for 
which  P-I  curves  might  be  formulated  was  initiated  by  reviewing  several 
different  Army  weapon  system  configurations.  Drawings  of  the  UH-1B 
helicopter,  the  S-280  electronic  equipment  shelter,  and  the  M-113  per¬ 
sonnel  carrier  were  studied  to  determine  what  simple  structural  elements 
comprising  these  systems  or  other  systems  might  be  considered  as  repre¬ 
sentative  or  typical  items  for  analysis.  During  the  reviewing  period, 
particular  attention  was  given  to  the  geometry  of  both  primary  and 
secondary  structure,  such  as  exterior  skin  panels  which  would  be  exposed 
to  blast  effects  and  interior  back-up  structure. 

Typical  panel  sizes,  expressed  as  length-to-width  ratios,  range 
from  2-4  for  the  M-113  personnel  carrier,  from  2.5-5  for  the  UH-1B 
helicopter,  and  as  much  as  2-8  for  the  S-280  electronic  equipment 
shelter.  Typical  panel  thicknesses  show  a  wide  variance  among  the 
three  weapon  systems  reviewed.  The  M-113  typical  panels  may  be  charac¬ 
terized  as  single-layered  thick  flat  plates  of  homogeneous  material, 
ranging  in  thickness  from  31.8  -  44.5  mm.  The  UH-1B  helicopter  panels, 
on  the  other  hand,  are  single-layered  thin  curved  plates  of  homogeneous 
material,  varying  in  thickness  from  .5  -  1.3  mm.  For  the  S-280  shelter, 
typical  skin  panels  are  of  heterogeneous  sandwich  construction  with  a 
typical  thickness  on  the  order  of  52  mm. 

The  primary  panel  support  or  back-up  structure  for  the  S-280 
shelter  is  a  brace  in  the  form  of  a  hat-section  extrusion.  Typical  of 
aircraft  semi-monocoque  construction,  the  UH-1B  back-up  structure,  con¬ 
sisting  of  frames,  longerons,  and  stringers,  is  comprised  of  a  variety 
of  extrusions  of  varying  length  and  cross-sectional  geometry.  Due  to 
the  thick  plate,  welded  construction  employed  in  the  M-113  personnel 
carrier,  secondary  or  back-up  structure  is  virtually  non-existent. 
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The  choice  of  the  range  of  geometrical  parameters  for  the  various 
types  of  structural  models  to  be  analyzed  was  relatively  uncomplicated 
except  for  the  beam  structural  models.  For  the  cylinder,  the  response 
is  a  function  of  two  geometrical  ratios,  the  length  to  radius  ratio  and 
the  radius  to  thickness  ratio.  For  flat  plates,  the  response  is  a  func¬ 
tion  of  the  length-to-width  and  width-to-thickness  ratios  of  the  plate. 

For  straight  beams,  however,  the  response  depends  upon  the  length  and 
the  cross-sectional  configuration. 

As  previously  mentioned,  the  backup  structures  for  the  Army  weapon 
systems  reviewed  consisted  of  a  variety  of  structural  elements  such  as 
braces,  stringers,  longerons,  etc.,  all  of  which  may  be  generally 
classified  as  structural  beams.  The  structural  beams  which  had  been 
looked  at  in  the  preliminary  review  task  had  cross  sections  such  as 
depicted  in  Figure  3.1  and  were  typical  of  standard  extruded  members. 
Generally  speaking,  all  of  these  extrusions  may  be  typified  as  either 
two-  or  three-area  cross  sections,  as  indicated.  Both  the  two-  and 
three-area  cross  sections  have  a  top  horizontal  area  Ax  attached  to  a 
vertical  area  Ax,  although  the  bottom  horizontal  area  A3  may  or  may  not 
exist  depending  upon  the  configuration  used.  Also,  as  indicated  in 
Figure  3.1,  either  area  A2  or  area  A3  may  be  comprised  of  more  than  one 
element.  The  vertical  distribution  of  the  individual  areas  of  a  beam 
cross  section  and  not  the  horizontal  distribution  of  the  areas  is  the 
important  factor  affecting  beam  vertical  bending  response. 

In  order  to  take  advantage  of  this  irrelevancy  of  the  horizontal 
distribution  of  individual  areas  of  beam  cross  sections,  it  was  decided 
to  synthesize  one  generic  beam  cross  section  which  could  be  used  to 
approximate  all  beam  cross  sections  of  interest.  The  synthesized  cross 
section  which  was  selected  for  the  straight  beam  model  P-I  analyses  is 
shown  in  Figure  3.2.  The  relative  depths  chosen  for  the  areas  Ax,  Ax, 
and  A3  were  based  on  a  review  of  the  pertinent  dimensions  of  typical 
standard  extrusions  tabulated  in  Ref.  4.  It  was  determined  from  the 
data  given  in  Ref.  4  that  a  typical  ratio  of  area  Ax  depth  to  either 
area  A^  or  area  A3  depth  was  on  the  order  of  18  to  1,  and  hence  this 
ratio  was  used  in  formulating  the  synthesized  cross  section.  Addition¬ 
ally,  it  was  determined  that  as  long  as  this  ratio  is  of  the  order 
indicated,  the  beam  response  will  be  insensitive  to  small  ratio  varia¬ 
tions. 

The  parameters  which  influence  the  beam  response,  exclusive  of  the 
parameters  defining  the  material  and  the  time-history  of  the  pressure 
loading,  are  the  two  area  ratios  A2/A^  and  A3/AX,  the  length-to-depth 
ratio  L/H,  the  value  of  H,  the  value  of  area  Ax,  and  the  loading  width 
parameter  VI.  Since  the  values  of  both  Ax  and  W  in  effect  singly  would 
alter  the  magnitude  of  the  beam  loading,  their  values  could  be  selected 
arbitrarily.  The  value  of  H  in  effect  also  alters  the  beam  loading  and, 
in  addition,  affects  the  time  scale.  For  input  simplification  purposes, 
therefore,  a  value  of  unity  was  chosen  for  W,  H,  and  wx  and,  as  can  be 
seen  from  Figure  3.2,  the  selection  of  a  value  for  wx  immediately  fixed 
the  value  of  Ax-  These  choices  also  result  in  Ax,  W,  and  H  (or  L)  being 
Involved  in  the  scaling  parameters  used  to  interpret  the  response 
results.  In  view  of  the  procedures  used  to  process  the  data  in  Section  4, 
this  scaling  is  implicit  rather  than  explicit. 

Ref.  4.  ALCOA  Structural  Handbook.  Aluminum  Company  of  America,  Pittsburgh, 
Pa.,  1958. 
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FIGURE  3.1  TYPICAL  EXTRUSION  CROSS  SECTIONS 


AREA  A|  0.05 Hw,  w( 


FIGURE  3.2  SYNTHESIZED  BEAM  CROSS  SECTION 


An  examination  of  the  extrusion  data  in  Ref.  4,  for  both  the  two- 
area  and  three-area  cross  sections,  was  then  made  to  determine  what 
ranges  of  the  area  ratios  indicated  on  Figure  3.2  should  be  considered. 
After  having  selected  specific  combinations  of  area  ratios  for  both  two- 
and  three-area  beam  cross  sections,  these  beam  cross-sectional  con¬ 
figurations  were  then  combined  with  selected  beam  lengths  to  define  the 
limits  of  the  beam  model  parameters  to  be  analyzed.  In  addition,  for 
curved  beams  various  typical  radii  of  curvature  were  selected  as  an 
additional  variable  parameter. 

3. 3  Model  Boundary  Conditions  and  Materials  of  Construction 

The  boundary  conditions  selected  for  the  structural  model  response 
analyses  were  considered  from  the  standpoint  of  actual  construction 
techniques,  practical  simulation  by  mathematical  formulation,  and  limi¬ 
tations  of  the  existing  response  codes.  Recognizing  the  fact  that  the 
actual  support  offered  to  individual  beams,  panels,  and  cylinders 
comprising  the  types  of  actual  structural  systems  considered  would 
probably  be  between  the  mathematical  representations  of  simple  supports 
and  clamped  supports,  it  was  adjudged  that  a  clamped  support  is  probably 
a  closer  approximation  to  reality  than  is  a  simple  support.  The  deci¬ 
sion  was  made,  therefore,  to  employ  a  rigidly-fixed  support  condition 
for  all  the  types  of  models  to  be  analyzed. 

All  of  the  previously  mentioned  weapon  systems  which  were  reviewed 
for  purposes  of  selecting  typical  structural  elements  for  blast  response 
analysis  utilize  different  alloys  of  aluminum  for  the  primary  and  sec¬ 
ondary  structural  elements.  The  M-113  personnel  carrier  plates  are 
formulated  from  5083  aluminum  alloy  material,  the  S-280  electronic 
equipment  shelter  employs  6061  aluminum  alloy  both  as  a  part  of  the 
outer  panel  sandwich  construction  and  for  the  back-up  structure  extru¬ 
sions,  and  the  UH-1B  helicopter  uses  2024  aluminum  alloy  for  the  back-up 
structure  and  both  2024  aluminum  alloy  and  AZ31B  magnesium  alloy  for  the 
outer  skin  panels.  Another  commonly  used  aluminum  alloy  for  which  good 
stress-strain  data  are  readily  available  is  7075  aluminum  alloy.  The 
decision  made  jointly  by  BRL  and  Raman  AviDyne  personnel  on  the  materials 
to  be  used  in  the  blast  response  analyses,  therefore,  was  to  consider 
one  2024  and  one  7075  aluminum  alloy.  Accordingly,  2024-T4  and  7075-T6 
were  selected  as  appropriate  aluminum  alloys  for  which  reasonable 
amounts  of  material  properties  data  were  readily  available. 

In  both  the  preliminary  response  analyses  and  the  final  vulnera¬ 
bility  code,  then,  structural  elements  comprised  of  either  2024-T4  or 
7075-T6  aluminum  alloy  and  rigidly  fixed  at  both  ends  have  been  assumed. 

3.4  Aerodynamic  Loading  Functions 

3.4.1  Loading  on  Circular  Cylinders 

Due  to  reflection  effects  which  occur  when  a  cylinder  is 
loaded  side-on  rather  than  end-on,  the  side-on  orientation  is  more 
critical  for  a  given  free-field  overpressure.  Hence,  the  side-on 
loading  condition  will  be  addressed  here. 
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The  side-on  loading  is  characterized  by  reflection  of  the 
shock  wave,  producing  high  pressures  which  are  subsequently  relieved  by 
diffraction  effects.  A  semi -empirical  model  for  the  loading  on  aircraft 
fuselages  is  given  in  Ref.  5.  Loading  on  a  cylinder  can  easily  be 
derived  from  this  semi-empirical  model.  It  turns  out  that  the  time 
constant  of  the  loading  through  the  reflection  and  diffraction  periods 
is  proportional  to  the  cylinder’s  radius.  The  time  constant  of  the 
structural  response  of  the  cylinder  is  also  proportional  to  the  radius. 
Hence,  if  response  only  during  the  reflection-diffraction  period  is 
considered,  the  physical  size  of  the  cylinder  makes  no  difference;  only 
ratios  are  required. 

The  semi -empirical  method  of  Ref.  5  remains  valid  only  as  long 
as  the  blast-wave  duration  is  long  in  comparison  with  the  reflection- 
diffraction  period.  There  appears  to  be  no  reliable  method  available 
for  blast-wave  durations  which  are  not  reasonably  long  in  comparison 
with  the  reflection-diffraction  period.  Furthermore,  structural  response 
runs  indicate  that  the  peak  response  of  a  cylinder  occurs  during  or 
fairly  shortly  after  the  reflection-diffraction  period.  Hence,  with  the 
restriction  that  the  blast-wave  duration  must  be  long  with  respect  to 
the  reflection-diffraction  period,  it  makes  little  difference,  as  far  as 
the  structural  response  is  concerned,  how  long  the  blast-wave  duration 
is.  With  little  error,  then,  the  blast  wave  may  be  taken  as  a  step 
function;  that  is,  with  infinite  positive  duration.  Thus,  the  present 
report  will  be  concerned  only  with  infinite-duration  blast  waves  for 
cylinders.  It  must  be  emphasized  that  the  development  and  verification 
of  a  loading  model  for  short  blast-wave  durations  would  require  reexami¬ 
nation  of  the  procedure.  In  the  meantime,  it  can  be  observed  that  use 
of  the  infinite-duration  blast  loading  will  result  in  overestimating 
the  damage,  but  probably  not  by  a  significant  amount  for  most  blast 
response  problems  of  practical  interest. 

Earl if  r,  it  was  indicated  that  the  structural  response  of  a 
cylinder  (in  terms  of  dimensionless  quantities,  such  as  strain)  is 
independent  of  the  size  of  the  cylinder  during  the  reflection-diffraction 
period.  The  restriction  to  infinite  duration  blast  waves  extends  this 
conclusion  to  the  entire  loading  period.  Hence,  the  computer  runs  can 
be  made  for  the  actual  cylinder  loading,  rather  than  for  an  arbitrary 
load-time  history  which  would  require  interpretation  through  pressure- 
impulse  considerations  for  the  actual  loading.  Stated  differently,  the 
present  study  of  cylinders  is  restricted  to  a  single  loading  time  his¬ 
tory;  hence,  no  requirement  exists  for  estimating  the  response  to  one 
loading  time  history  from  calculated  response  to  a  different  time 
history  -  the  actual  time  history  can  be  used  to  calculate  the  response 
in  the  first  place. 


Ref.  5.  NOVA  -  A  Digital  Computer  Program  for  Calculating  the  Response 
of  Aircraft  to  the  Overpressure  from  a  Nuclear  Weapon.  Vol.  I  - 
Theory  and  Program  Description.  Prepared  by  Katnan  AviDyne  for 
'*  th  AFWL,  AFWL-TR-72-115,  Vol.  I,  July  1973. 
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Attention  is  accordingly  directed  toward  defining  the  loading 
time  history  for  a  blast  wave  of  infinite  duration  based  on  the  semi- 
empirical  method  for  fuselage  loadings  in  Ref.  5.  The  method  of  Ref.  5 
is  more  comprehensive  than  required  here,  since  it  treats  arbitrary 
blast-structure  orientations  rather  than  just  side-on  and  because  it 
deals  with  non-cylindrical  shapes  which  may  have  a  forward  velocity. 

The  forward  velocity  makes  no  difference  for  a  right  circular  cylinder. 

The  typical  time  history  of  the  pressure  on  the  cylinder  is 
shown  in  Figure  3.3.  The  shock  front  arrives  at  time  ta,  which  is  a 
function  of  the  angle  0,  defined  in  the  insert  in  Figure  3.3.  With  zero 
time  defined  as  the  time  at  which  the  shock  front  first  reaches  any 
point  on  the  cylinder, 


R(l-cos9) 
a  V 

s 


(3.1) 


where 


R  is  the  cylinder  radius 
Vg  is  the  shock  front  velocity 


The  blast  wave  may  be  characterized  by  the  overpressure  associated  with 
it,  Apg.  The  shock  strength  is  defined  as 
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where 


Pq  is  the  ambient  pressure 

The  shock  velocity  may  be  found  from  the  shock  strength  using  the 
Rankine-Hugoniot  relation 
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where  cq  is  the  ambient  speed  of  sound. 
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The  time  t<j  in  Figure  3.3  is  the  time  at  which  the  reflection- 
diffraction  effects  end. 

From  Ref.  5 


The  loading  is  taken  as  linear  from  t  to  t  . 

a  d 
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where  a  is  the  speed  of  sound  behind  the  shock  front. 


and  is  given  by 
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Equation  (3.4)  results  in  t<j  =  0  for  0=0.  This  is  not  correct,  and  to 
avoid  this  situation.  Ref.  5  introduces  a  minimum  reflection-diffraction 
period  given  by 
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Having  defined  the  times  associated  with  the  cylinder  loading  time 
history,  the  only  remaining  variables  to  be  defined  are  the  pressures  at 
times  tfl  and  t,,  the  pressure  being  constant  for  t^tj  due  to  the  assump¬ 
tion  of  a  blast  wave  of  infinite  duration.  From  Ref.  5,  the  loading  at 

time  t  is 
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where  r  is  the  reflection  factor  for  normal  reflection  of  a  shock  wave, 
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Ref.  5  as 


Finally,  the  pressure  loading  at  time  t  , 

-  d 


may  be  found  from 


p(td,0)  =  Apg  +  Cpq 


(3.9) 


where 

Cp  is  a  pressure  coefficient 

q  is  the  dynamic  pressure  behind  the  shock  front 


The  pressure  coefficient  is  represented  in  Ref.  5  as 
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The  dynamic  pressure  may  be  expressed  in  terms  of  the  shock  strength,  by 
using  the  Rankine-Hugoniot  relations,  as 
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where  p  is  the  ambient  density, 
o 

This  completes  the  definition  of  the  pressure  loading  on  a  cylinder  as 
utilized  in  both  the  structural  model  response  analysis  and  the  vulnera¬ 
bility  code. 

3.4.2  Loading  for  Beams  and  Plates 

The  considerations  which  led  to  the  use  of  the  actual  side-on 
loading  for  cylinders  are  not  valid  when  beams  and  plates  are  addressed. 
For  example,  the  duration  of  the  reflection-diffraction  period  is  a 
function  of  the  size  of  the  structure  of  which  the  beam  is  a  pant,  and 
this  size  is  not  defined  by  the  size  of  the  beam.  The  problem  of  defin¬ 
ing  the  pressure  loading  for  blast-wave  durations  which  are  not  long  in 
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comparison  with  the  reflection-diffraction  period  remains;  however,  the 
beam  may  have  a  long  enough  response  time  constant  so  that  the  peak 
response  may  occur  well  after  the  reflection-diffraction  period.  Hence, 
the  assumption  of  an  infinite  duration  blast  wave  is  not  viable  for 
beams  and  plates. 

Accordingly,  the  loading  time  history  used  for  the  structural 
model  response  analyses  to  establish  the  pressure-impulse  curves  was  a 
simple  exponential  form, 


p(t) 


(3.12) 


where  p^  and  tG  are  pressure  and  time  parameters  which  permit  different 
values  of  pressure  and  impulse  to  be  examined.  The  time  history  given 
by  Eq.  (3.12)  is  similar  to  that  which  will  actually  be  encountered,  in 
that  it  jumps  to  its  peak  at  time  zero  and  subsequently  decays  in  an 
exponential  fashion.  For  the  pressure-impulse  method  to  work,  it  is 
important  that  the  loading  time  history  on  which  the  P-I  curves  are 
based  be  similar  to  the  loading  time  histories  for  which  the  P-I  curves 
will  be  used. 

For  the  vulnerability  code,  the  shock  front  is  assumed  to 
encounter  the  structure  side  on;  that  is,  with  the  face  of  the  structure 
containing  the  beam  or  panel  oriented  parallel  to  the  shock  front.  This 
is  at  least  close  to  being  the  most  critical  orientation.  The  vulner¬ 
ability  code  then  has  imbedded  in  it  a  model  for  the  side-on  blast- 
wave  loading  of  a  flat  (or  nearly  flat)  structure.  This  model  is  taken 
from  Ref.  6. 


The  pressure  loading  at  time  zero  is  the  same  as  that  on  a 
cylinder  at  0=0;  that  is,  from  Eq.  (3.7) 


p(0>  =  Ap  r  (3.13) 

F 

where  r  is  given  by  Eq.  (3.8),  To  define  the  reflection-diffraction 
period, ptwo  distances,  SS  and  S  are  defined  in  Ref.  6.  The  first,  SS, 
is  the  distance  from  the  point  in  question  (taken  here  to  be  the  center 
of  the  beam  or  plate)  to  the  nearest  free  edge.  A  free  edge  as  used 
here  is  represented  by  a  position  on  the  structure  at  which  the  struc¬ 
ture  is  perpendicular  to  the  shock  wave.  For  example,  in  the  case  in 
which  the  front  face  of  a  box  structure  on  the  ground  is  being  examined, 
the  intersections  of  the  roof  and  side  walls  with  the  front  face  would 
be  free  edges.  The  distance  SS  is  used  to  define  the  time  at  which 
diffraction  effects  begin,  t 


Ref.  6.  Dunn,  B.J.,  Lavagnino,  A.,  and  White,  J.J.,  Parametric  Analysis 
of  Nuclear  Weapon  Air  Blast  Effects  on  a  Model  Target,  BDM/W- 
74-009TR,  Braddock,  Dunn  and  McDonald,  Inc.,  December  1973. 
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(3.14) 


s 

where  V  ,  the  shock  velocity,  is  given  by  Eq.  (3.3). 

The  second  distance  defined  in  Ref.  6,  S,  is  the  longest 

distance  from  any  point  on  the  face  of  the  structure  being  considered 

to  the  farthest  free  edge.  The  distance  S  is  used  to  define  the  time  at 

which  diffraction  effects  end,  t  . 

s 


t 


s 


(3.15) 


From  time  zero  to  time  tss,  the  loading  will  be  taken  to  be 
linear.  In  accordance  with  Ref.  6,  the  loading  at  time  t  is 


p(tss> 


r  Ap (t  ) 
p  ss 


(3.16) 


where  Ap(t)  is  the  overpressure  associated  with  the  blast  wave  at  time 
t.  It  is  seen  that  the  decay  in  loading  from  time  zero  to  time  tss 
is  due  to  blast  wave  decay  rather  than  to  diffraction  effects. 

From  time  tgg  to  time  tg,  diffraction  effects  cause  the 
loading  to  diminish  to  the  steady  state  pressure  loading  at  time  t  , 
given  by 

p(tg)  =  Ap  (tg)  +  q(tg)  (3.17) 


where  q(t)  is  the  dynamic  pressure  associated  with  the  blast  wave  at 
time  t  and  a  drag  coefficient  of  unity  has  been  assumed.  From  time  tgg 
to  time  t  ,  the  pressure  loading  is  taken  to  be  linear.  The  form  of 
Eq.  (3,17)  is  applicable  for  any  time  subsequent  to  tg;  that  is. 


p(t)  =  Ap (t)  +  q(t),  t_>tg 


(3.18) 


Equation  (3.11)  expresses  the  dynamic  pressure  in  terms  of  the  shock 
strength,  which  is  defined  by  Eq.  (3.2).  The  relationship  thus  expressed 
between  dynamic  pressure  and  overpressure  is  properly  valid  only  at  the 
shock  front;  it  may  be  used  as  an  approximation  at  points  behind  the 
shock  front;  hence, 
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(3.19) 


Substituting  Eq.  (3.19)  into  Eq .  (3.18), 


p(t)  =  Ap(t) 


25 


Ap(t) 


1  + 
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~  P  c 
2  o  o 


7,7  +  AElOv 

Po  I 


(3.20) 


This  completes  the  development  of  the  pressure  loading  models  for  plates 
and  beams  in  both  the  structural  model  response  analyses  and  the  vulner¬ 
ability  code. 

3 . 5  Response  Damage  Criteria 

The  structural  model  response  or  damage  levels  which  were  used  as 
response  criteria  were  the  desired  iso-damage  levels  selected  by  BRL  for 
which  P-1  curves  would  be  formulated.  Three  distinct  damage  levels  weri 
specified  and  these  levels  correspond  to  light,  moderate,  or  heavy 
damage,  respectively. 

Since  the  response  analyses  required  a  more  precise  definition  of 
response  criteria  than  light,  moderate,  or  heavy  damage,  the  light 
damage  level  was  specified  as  inelastic  response  strain  of  3%;  the  heavy 
damage  level  was  specified  as  inelastic  response  strain  corresponding  to 
material  rupture;  and  the  moderate  damage  level  was  specified  as  u 
response  corresponding  co  an  inelastic  strain  level  halfway  between 
light  and  heavy  damage. 

3.6  Methods  of  Response  Analysis 

Having  described  the  types  of  structural  elements  requiring  analy¬ 
sis,  the  model  boundary  conditions  and  materials  of  construction,  the 
aerodynamic  loading  functions,  and  the  response  damage  criteria,  atten¬ 
tion  may  now  be  given  to  a  description  of  the  methods  of  response 
analysis . 


3.6.1  Cylinder  Response  Analysis 

As  previously  indicated  in  Subsection  3.4,  the  cylinder  response 
investigations  were  limited  to  infinite  duration  blast-wave  loadings.  Hence, 
no  P-I  data  points  per  se  were  obtained  in  the  cylinder  response  analyses, 
only  pressures  defined  by  the  previously  stated  forcing  function  and  cor¬ 
responding  to  the  three  levels  of  damage  criteria.  The  method  used  to  ana¬ 
lyte  the  response  of  the  cylinder  models  was  the  DEPICS  computer  code  (Ref.  7) 
which  was  developed  to  predict  the  dynamic  response  of  multi-layered  cylin¬ 
drical  or  conical  bodies  undergoing  large  non-linear  elastic  or  inelastic 
deformations. 


Ref.  7.  Mente,  L.J.,  The  Dynamic,  Elastic-Plastic,  Large  Displacement  Re¬ 
sponse  of  Buckling  Sensitive  Cylindrical  Shells  to  Blast-Type 
Loadings,  —  Part  I  Analytical  Formulation.  AMC-2-68-(T) ,  Kaman 
AviDyne  TR-53,  Kaman  Sciences  Corporation,  August  1968. 
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It  had  originally  been  anticipated  that  DEP1CS  would  be  run  in 
an  iterative  mode  to  permit  the  program  to  seek,  the  pressure  which  would 
result  in  a  given  level  of  response  for  a  specific  combination  of  the 
cylinder  parameters.  However,  since  DEPICS  would  have  required  a  con¬ 
siderable  amount  of  computer  time  for  each  iterative  run  desired,  it  was 
subsequently  decided  chat  it  would  be  simpler  to  run  DEPICS  in  a  straight 
response  mode  to  obtain  a  number  of  strain  response  values  corresponding 
to  selective  input  values  of  the  forcing  function  pressure.  By  follow¬ 
ing  this  approach  for  each  combination  of  cylinder  parameters,  it  was 
possible  to  construct  a  number  of  curves  relating  cylinder  response 
strain  to  forcing  function  pressure.  Each  individual  curve  was  then 
entered  consecutively  with  the  three  specified  levels  of  inelastic 
strain  criteria  and  the  corresponding  pressure  read  off. 

The  cylinder  model  configurations  analyzed  included  length-to- 
diameter  ratios  between  1  and  3  and  radius-to-thickness  ratios  between 
50  and  300.  These  limits  are  not  expected  to  be  exceeded  in  actual 
weapon  system  configurations. 

3.6.2  Beam  Response  Analyses 

The  beam  response  investigations  were  considerably  more  com¬ 
plicated  than  the  cylinder  analyses  due  to  the  larger  number  of  param¬ 
eters  involved  and,  additionally,  due  to  the  fact  tha'  P-I  data  were 
being  sought. 

The  method  of  analysis  used  to  obtain  beam  response  was  a 
modified  version  of  the  beam  option  of  Kaman  AviDyne's  blast  response 
code,  NOVA  (Ref.  5).  The  modifications  consisted  of  two  types;  first, 
basic  modifications  to  the  program  to  permit  accurate  calculation  of 
stresses  and  strains  near  the  clamped  edge,  and  second,  modifications 
designed  to  streamline  the  calculation  of  P-I  curves. 

Consid*  ring  the  first  type  of  modif ication,  as  indicated  in 
Subsection  3.3,  clamped  edge  conditions  were  assumed.  Comparison  of 
results  from  the  beam  code  with  theoretical  static  elastic  analyses 
indicated  that  stresses  near  the  clamped  edge  were  badly  underpredicted, 
particularly  for  relatively  thin  beams.  The  problem  arises  basically 
due  to  the  fact  that  a  very  high  stress  gradient  occurs  near  the  clamped 
edge,  and  the  beam  program  gave  insufficient  definition  to  follow  this 
gradient.  The  beam  program  breaks  the  beam  down  into  weightless  beam 
segments  connecting  discrete  mass  points.  The  spacing  is  usually  taken 
to  be  fairly  even  to  conserve  computer  time,  since  the  integration  time 
interval  which  can  be  used  is  proportional  to  the  length  of  the  shortest 
segment.  It  was  found  that  the  addition  of  many  closely-spaced  mass 
points  near  the  clamped  edge  permitted  accurate  calculation  of  the 
static  edge  stresses.  This  procedure  could  not  be  used  for  dynamic 
analyses,  however,  because  the  aforementioned  dependence  of  the  inte¬ 
gration  time  interval  upon  segment  length  would  have  resulted  in 
prohibitively  long  computer  running  times.  The  solution  to  this  problem 
involves  extensive  modification  to  the  computer  code  so  that  the  beam 
segment  near  the  clamp  is  handled  in  a  pseudo-static  manner.  This 
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approach,  which  is  described  in  detail  in  Ref.  8,  resulted  in  computer 
running  times  approximately  two  orders  of  magnitude  less  than  would  have 
been  required  by  adding  sufficient  mass  points  while  predicting  vir¬ 
tually  the  same  response. 

The  second  type  of  modification  to  NOVA  was  aimed  at  streamlining 
the  calculation  of  P-I  curves.  NOVA  was  originally  formulated  to  cal¬ 
culate  the  overpressure  response  corresponding  to  a  preselected  range 
and  weapon  yield,  and  the  overpressure  time-history  corresponding  to  the 
input  range  and  yield  was  internally  calculated  within  NOVA  through  the 
employment  of  a  self-contained  blast  model.  In  order  to  obtain  more 
definitive  control  over  the  magnitude  and  duration  of  the  forcing 
function,  NOVA  was  modified  for  the  P-I  analysis  problem  to  bypass  the 
original  self-contained  blast  model,  and  the  basic  beam  and  plate 
forcing  function  time-history  expressed  by  Eq.  (3.12)  was  inserted  into 
NOVA.  Additionally,  the  program  was  modified  to  iterate  on  the  magni¬ 
tude  of  the  forcing  function  to  obtain  a  preselected  response  criterion 
corresponding  to  a  preselected  value  of  the  forcing  function  duration. 

The  previously  described  modifications  to  the  original  NOVA  computer 
code  ultimately  resulted  in  a  revised  computer  code  called  KADBOP,  an 
acronym  for  Kaman  AviDyne  Dynamic  response  of  Beams  or  Plates,  which  was 
specifically  tailored  to  the  requirements  of  the  planned  P-I  investigations. 

In  order  to  assist  the  reader  in  understanding  how  the  basic 
beam  option  of  KADBOP  was  intended  to  be  used,  a  brief  explanation  of 
the  procedures  to  be  followed  in  obtaining  one  beam  P-I  data  point  will 
be  presented.  First,  for  the  particular  beam  material  selected,  a 
particular  value  of  L/H  is  chosen.  Since,  as  previously  indicated,  the 
input  value  of  H  was  selected  as  unity  for  simplification  purposes,  the 
length  of  the  beam  between  clamped  supports  is  established  and  is  then 
appropriately  modeled  to  meet  the  KADBOP  beam  response  input  require¬ 
ments.  Next,  values  for  the  area  ratios  shown  in  Figure  3.2  are  selec¬ 
ted  and  the  input  values  W£  and  W3  calculated.  The  material  stress- 
strain  data  and  the  desired  three  levels  of  response  criteria  are  then 
specified  in  the  program  input.  Finally,  a  value  of  tQ  and  an  initial 
value  of  pM  for  the  forcing  function  are  chosen  and  included  in  the 
input  set.  At  this  point,  the  KADBOP  beam  response  program  has  all  the 
data  required  and  is  ready  to  calculate  beam  dynamic  response. 

Since  each  point  on  a  P-I  curve  defines  a  pressure  and  impulse 
combination  corresponding  to  a  specific  level  of  damage  or  response,  the 
program  has  to  be  run  In  the  iterative  option  to  determine  the  required 
pressure  loads.  Basically,  in  the  iterative  option  of  KADBOP,  the 
program  holds  the  input  loading  duration  parameter,  tc,  constant  and 
adjusts  the  input  loading  force  parameter,  p  ,  until  the  response 
matches  the  input  response  criterion.  Therefore,  using  the  input  value 


Ref.  8.  NOVA-2  -  A  Digital  Computer  Program  for  Analyzing  Nuclear 

Overpressure  Effects  on  Aircraft,  Prepared  by  Kaman  AviDyne 
for  the  AFWL,  KA  TR-128,  (to  be  published). 
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of  tQ  and  the  three  levels  of  strain  criteria,  the  program  computes 
three  force  factors  which,  when  multiplied  by  the  input  </alue  of  pQ, 
determines  one  pressure  data  point  for  each  of  the  three  P-I  damage 
curves  being  constructed.  Once  the  pressure  data  point  is  calculated, 
it  is  a  simple  matter  to  determine  the  corresponding  impulse,  since 


I 


p(t)dt. 


(3.21) 


This  same  procedure  would  normally  be  repeated  for  other  input  values  of 
tQ  until  enough  P-I  data  points  are  available  to  construct  the  desired 
P-I  iso-damage  curves. 

The  original  beam  option  of  KADBOP  just  described  was  intended 
to  be  applied  to  the  various  combinations  of  beam  geometry  parameters 
which  had  been  selected  as  typical  of  the  beam  geometries  which  would 
probably  be  encountered  in  real  weapon  system  configurations.  The 
ranges  of  beam  parameters  selected  for  analysis  were:  beam  length-to- 
depth  ratios  of  15  to  75;  A2/A^  area  ratios  of  .3  to  3.6;  and  A^/Aj 
area  ratios  of  0  to  1.5. 

It  was  further  intended  that  the  P-I  data  points  emerging  from 
the  KADBOP  response  analysis  would  be  used  to  construct  P-I  iso-damage 
curves  for  various  combinations  of  the  specified  parameters  and  that 
these  curves  would  then  be  imbedded  in  the  proposed  vulnerability  code 
as  tables  of  discrete  values.  The  use  of  the  proposed  vulnerability 
code  for  geometry  configurations  intermediate  to  the  imbedded  geometries 
would  then  have  required  multiple  interpolations.  It  will  subsequently 
be  indicated  in  Section  4  that  the  proposed  interpolation  posed  a  prob¬ 
lem  and  that  a  more  accurate  method  was  adopted  for  imbedding  P-I  data 
in  the  vulnerability  code. 

As  can  be  seen  from  the  large  number  of  parameters  selected 
and  the  further  consideration  of  the  additional  parameters  required  for 
curved  beams,  the  number  of  P-I  data  points  to  be  calculated  by  the 
Iterative  procedure  previously  described  was  quite  extensive  when  all 
the  possible  combinations  of  the  beam  parameters  were  considered.  Con¬ 
sequently,  as  will  be  described  in  the  next  subsection,  a  method  was 
sought  which  would  reduce  the  number  of  computer  runs  required  to  accu¬ 
mulate  this  large  mass  of  P-I  data.  Before  considering  this  aspect  of 
the  P-I  data  development,  however,  the  method  of  response  analysis  for 
plates  should  be  discussed. 

3.6.3  Plate  Response  Analysis 

As  indicated  above,  extensive  modification  of  an  existing  beam 
program  was  necessary  in  order  to  permit  accurate  calculation  of  stresses 
and  strains  near  the  clamped  boundary.  Similarly,  it  was  found  that 
plate  codes  which  were  available  could  not  calculate  accurate  edge 


stresses  near  a  clamped  boundary  in  a  practical  manner.  (A  very  recent 
code,  reported  in  Ref.  8,  does  have  this  capability,  but  was  not  avail¬ 
able  in  time  to  be  used  for  the  present  program.) 

One  alternative  method  available  was  to  use  the  beam  code, 
which  would  give  essentially  the  correct  result  for  moderate  to  high 
aspect  ratio  plates,  and  to  modify  the  beam  results  for  application  to 
squarish  plates  using  results  of  theoretical  static  elastic  analyses  of 
plates.  The  other  alternative  was  to  use  an  available  plate  code  and  to 
recognize  that  the  edge  strains  would  not  be  accurately  calculated.  The 
decision  was  made  in  favor  of  the  first  alternative,  primarily  because 
it  seemed  very  desirable  to  maintain  consistency  between  the  beam  and 
plate  results. 

3.6.4  Comments  on  Analyses 

As  stated  earlier,  the  analyses  are  all  based  on  clamped 
boundary  conditions.  Although  this  boundary  condition  is  probably  the 
most  appropriate  one  to  use,  the  error  introduced  by  replacing  the 
actual  boundary  condition,  whatever  it  is,  by  an  ideal  clamp  is  not 
known. 


Further,  the  use  of  a  clamped  boundary  for  beams  forces  the 
analyst  to  consider  the  difficult  problems  associated  with  the  very  high 
strain  gradient  near  the  clamped  boundary.  Questions  arise  as  to  the 
validity  of  the  analysis  near  the  boundary  when  applied  to  an  actual 
structure.  The  actual  structure  may  be  stronger  than  the  idealized 
structure  which  was  analyzed  because  the  actual  edge  condition  is  not  an 
ideal  clamp.  On  the  other  hand,  it  may  be  weaker  due  to  the  presence  of 
edge  fastenings  which  result  in  local  stress  concentrations.  Finally, 
there  is  the  question  of  how  to  apply  tensile-test  rupture  data  to  a 
case  in  which  the  high  strain  is  very  localized  and  results  largely  from 
bending  rather  than  axial  load. 

The  analysis  methods  employed  are  geometrically  as  well  as 
materially  non-linear;  they  are  capable  of  predicting  overall  buckling 
of  a  curved  beam,  for  example.  They  will  not  predict  local  crippling, 
such  as  might  occur  in  an  I-beam  flange.  Thus,  crippling  is  ignored  in 
the  present  work.  Since  crippling  will  weaken  a  beam,  neglecting 
crippling  will  lead  to  an  underprediction  of  damage. 

3.7  P-I  Curve  Hyperbolic  Approach 

The  previous  subsection  described  the  methods  of  analysis  and  com¬ 
puter  codes  which  were  employed  to  initiate  the  studies  of  the  non¬ 
linear  dynamic  response  of  beams,  plates,  and  cylinders  in  order  to 
generate  P-I  data.  As  previously  indicated,  the  cylinder  analyses 
yielded  curves  of  pressure  versus  response  strain  for  various  cylinder 
configuration  parameters  and,  as  such,  did  not  require  formulation  of 
the  usual  P-I  iso-damage  curves.  The  beam  and  panel-beam  analyses  as 
described,  however,  do  yield  sets  of  discrete  P-I  data  points  for  what¬ 
ever  combination  of  beam  configuration  parameters  are  being  considered. 
Figure  3.4  portrays  a  typical  light  damage  curve  for  one  beam  material, 
one  specific  L/H  value,  one  specific  value  of  A2/A^,  and  one  value  of 
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FIGURE  3.4  TYPICAL  BEAM  P-I  CURVE 


A^/A^.  The  construction  of  this  curve  represents  a  considerable  number 
of  discrete  P-I  data  points  and,  consequently,  required  a  considerable 
number  of  computer  response  runs.  It  rapidly  became  evident  during  the 
course  of  the  beam  and  plate  studies  that  the  total  number  of  required 
computer  response  runs  and  the  attendant  computer  running  time  would 
reach  staggering  proportions  in  order  to  fulfill  all  the  planned  para¬ 
metric  investigations.  Consequently,  a  method  was  searched  for  which 
would  provide  the  required  P-I  data  but  with  much  less  effort. 

It  has  been  previously  recognized  by  other  investigators  that  many 
P-I  curves  are  hyperbolic  in  nature  and  approach  asymptotic  values  of 
both  pressure  and  impulse.  On  the  assumption  that  a  given  P-I  iso¬ 
damage  curve  is  a  true  hyperbola,  it  is  possible  then  to  take  advantage 
of  the  known  equation  for  a  true  hyperbola  to  derive  several  useful 
mathematical  expressions. 

The  equation  for  a  hyperbola  is 


{k  - l)  (t  -  ^  ■ J 

where  PG  and  IQ  are  the  pressure  and  impulse  asymptotes  and  J  is  a 
constant.  The  asymptotes  are  found  by  KADBOP  by  using  a  step  function 
(tQ->-  ®  in  Eq.  (3.12))  for  the  pressure  asymptote  and  by  using  an  impulse 
(tQ  +0,  pm  ^  00  in  Eq.  (3.12))  for  the  impulse  asymptote.  Actually,  for 
the  impulse  asymptote,  the  initial  velocity  imparted  to  the  mass  points 
by  the  impulse  is  found  and  the  problem  is  then  simply  an  initial  con¬ 
dition  problem. 

In  order  to  determine  the  constant  J,  one  more  piece  of 
information  is  required.  The  hyperbola  given  by  Eq.  (3.22)  is  plotted 
in  Figure  3.5.  The  form  is  seen  to  be  symmetric  with  the  same  scale 
used  for  both  axes.  The  midpoint  or  knee  point,  occurs  at  an  angle  of 
45°  from  either  axis,  and  is  designated  as  point  K.  If  the  pressure  and 
impulse  at  point  K,  Pg  and  Ik,  are  defined,  substitution  of  these  values 
into  Eq.  (3.22)  will  allow  J  to  be  calculated.  Along  a  45°  line  on 
P  I 

Figure  3.5,  —  =  — •  By  integrating  Eq.  (3.12),  the  impulse  is  found 
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Hence,  the  equation  which  defines  the  45°  line  on  Figure  3.5  is 
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FIGURE  3.5  HYPERBOLIC  REPRESENTATION  OF  P-I  CURVE 
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Thus,  when  PQ  and  IQ  have  been  found,  tQ  can  be  calculated.  Using 
t^45)  an(j  an  arbitary  value  of  p  ,  the  response  level  which  occurs  at  a 
point  along  the  45°  line  may  be  determined.  By  varying  the  value  of  pM, 
the  point  at  which  the  desired  value  of  response  occurs  may  be  found, 
and  the  corresponding  value  of  p^  is  the  value  sought,  P^.  The  con¬ 
stant,  J,  may  then  be  calculated  as  indicated. 

As  a  check  on  the  accuracy  of  the  computed  hyperbola,  the  modified 
KADBOP  program  further  proceeds  to  calculate  the  response  of  one  pair  of 
P-I  values  on  either  side  of  the  knee  of  the  curve  and  compares  these 
values  with  the  strain  criterion  associated  with  the  particular  iso-damage 
curve  being  sought.  These  checks  demonstrated  the  validity  of  the 
assumption  of  a  hyperbolic  shape  for  the  P-I  curve. 

3.8  P-I  Curve  Usage 

Suppose  that  at  this  point  the  analyst  has  available  to  him  iso¬ 
damage  curves  for  light,  moderate,  and  heavy  damage,  and  that  one  set  of 
curves  is  for  beam  parameters  which  match  a  beam  he  wishes  to  analyze. 

The  analyst  may  want  to  determine  the  level  of  response  to  some  pressure 
time  history,  either  analytically  calculated  or  experimentally  measured. 

He  can  define  the  peak  pressure,  P,  and,  by  integration  over  the  posi¬ 
tive  duration,  the  impulse,  I,  and  can  thus  define  a  point  on  the  P-I 
plane.  By  noting  where  this  point  falls  relative  to  the  iso-damage 
curves  for  light,  moderate,  and  heavy  damage,  he  can  qualitatively 
define  the  desired  level  of  response;  for  example,  he  might  find  that 
the  response  will  be  between  light  and  moderate  damage. 

In  Section  2,  an  alternative  P-I  approach  based  on  work  by  Young- 
dahl  in  Ref.  2  and  modified  by  Cummings  and  Schumacher  in  Ref.  3  is 
mentioned.  Use  of  this  approach  is  a  bit  more  involved  than  the  straight 
P-I  approach  described  above.  The  same  iso-damage  curves  used  in  the 
above  example  will  be  used  here,  and  the  analyst  is  again  assumed  to 
have  a  pressure  time  history.  To  implement  the  modified  Youngdahl 
approach,  the  analyst  must  also  select  a  critical  pressure  value,  Pcr. 
Permissible  values  for  Pcr  in  the  present  application  range  from  zero  to 
the  value  of  the  pressure  asymptote  for  light  damage.  In  principal,  the 
value  of  PCr  should  be  related  to  the  value  of  static  pressure  which 
will  just  produce  inelastic  deformation  of  the  structure. 
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The  analyst  defines  an  effective  impulse  as 


t 

cr 

I  e  m  [  P<Odt  (3.24) 

o 


where  t  is  the  time  at  which  p(t)  decreases  to  Pcr.  It  is  assumed 
that  p(t)  is  of  the  form  considered  in  this  report;  that  is,  starting  at 
its  maximum  value  and  subsequently  decaying  in  a  more-or-less  expo¬ 
nential  manner.  The  effective  pressure  is  then  defined  as 


P 

e 


cr 

/ 


tp(t)dt 


(3.25) 


The  values  of  Pe  and  le  defined  using  Eqs.  (3.24)  and  (3.25)  cannot 
be  used  with  the  P-1  iso-damage  curves  directly.  Rather,  the  following 
procedure  must  be  used.  The  P-I  iso-damage  curves  are  based  on  the 
exponential  loading  given  by  Eq.  (3.12),  which  involves  the  two  param¬ 
eters  pm  and  tQ.  It  is  possible  to  determine  the  values  of  p^  and  tD 
which  would  yield  the  defined  values  of  Pe  and  Ie  for  the  given  value  of 
Pcr.  Having  pM  and  tQ,  the  pressure,  P,  is  immediately  available  as  pj^ 
and  the  impulse,  I,  is  available  through  integration  as  given  by  Eq.  (3.23). 
These  values  define  a  point  on  the  P-I  plane  which  may  be  used  as  described 
previously. 
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SECTION  4 

PROCESSING  DATA  FOR  VULNERABILITY  CODE 


4.1  Introduction 

The  previous  sections  have  described  the  codes  used  to  obtain  the 
non-linear  response  of  beams,  panels,  and  cylinders.  The  assumption 
that  the  P-I  curves  are  hyperbolas  has  been  introduced  and  justified. 

At  this  point,  then,  the  three  data  points  which  define  a  hyperbola  are 
available  for  a  number  of  beams,  panels,  and  cylinders  for  three  levels 
of  response.  The  next  task  is  to  process  these  data  for  use  in  a  vul¬ 
nerability  program. 

For  the  cylinders,  no  problem  exists.  The  number  of  parameters 
involved  is  so  small  (R/t  and  L/R)  that  numerical  interpolation  may  be 
conveniently  and  easily  used  to  obtain  desired  data  for  any  combination 
of  R/t  and  L/R  within  the  range  considered. 

The  situation  is  far  different  for  the  structures  analyzed  within 
the  beam  code.  Here  there  are  more  parameters  involved  (L/H ,  A^/A^, 

A /An  and  L/R),  with  the  consequence  that  insufficient  data  are  avail¬ 
able  to  cover  the  parameter  ranges  considered  in  the  detail  required  to 
permit  reasonably  accurate  numerical  interpolation.  Hence,  some  means 
of  fitting  the  data  was  sought  which  would  be  valid  over  the  entire 
region. 

It  would  be  possible  to  establish  fitting  functions  in  a  purely 
mathematical  fashion  utilizing  only  the  available  data.  However,  it  was 
deemed  preferable  to  draw  upon  knowledge  of  structural  behavior  to 
assist  in  arriving  at  appropriate  fitting  functions.  The  postulated 
procedure  was  to  utilize  static  elastic  solutions  to  guide  the  selection 
of  fitting  functions.  This  procedure  was  tried  first  for  straight 
beams,  with  excellent  results.  The  curved  beam  posed  a  much  more 
challenging  problem,  due  primarily  to  the  snap-through  phenomenon.  For 
this  reason,  the  fitting  functions  which  were  selected  fit  the  straight 
beam  data  considerably  better  than  the  curved  beam  data. 

The  static  elastic  analysis  is  described  in  Subsection  4.2  and  the 
reduction  of  the  results  to  the  straight  beam  case  is  given  in  Sub¬ 
section  4.3.  The  development  of  the  fitting  functions  is  presented  in 
Subsections  4.4  and  4.5  for  straight  and  curved  beams,  respectively. 

The  application  of  the  results  for  2024  A1  beams  to  7075  A1  and  to 
plates  is  described  in  Subsection  4.6,  and  the  accuracy  of  the  fitting 
functions  is  estimated  in  Subsection  4.7. 
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4. 2  Static  Elastic  Analysis 

The  static  elastic  problem  which  is  addressed  is  depicted  in  Figure 
4.1.  The  curvature  is  assumed  to  be  small,  with  the  result  that  the 
axial  reaction  force,  P,  is  equal  to  the  compressive  load  throughout  the 
length  of  the  beam.  The  equilibrium  equation  is 


d  (y-y  ) 

El  - - -  =  -Py 

dx 


(4.1) 


where 

El  is  the  bending  stiffness 

yQ  defines  the  initial  shape  of  the  beam 

and  the  remaining  variables  are  defined  by  Figure  4.1.  The  initial 
shape  is  assumed  to  be  parabolic,  which  is  somewhat  more  tractable  than 
a  circula-  shape,  leading  to 


The  solution  of  Eq.  (4.1)  is 


-2  1 

'wL  L  ; 

cos4>x  -  cos4>  1-x 

2P  "  2R^ 

<(>sin<J>  2 

where 


(4.3) 


x 


The  .;olution  given  by  Eq.  (4.3)  is  not  complete  because  the  value 
of  P  is  not  known.  The  deformation  shape  defined  by  Eq.  (4.3)  may  be 
used  to  determine  the  change  in  length  of  the  beam,  AL,  and  hence  the 
average  axial  strain.  The  relationship 


P  = 


(4.5) 


where  A  is  the  cross-sectional  area  of  the  beam  becomes,  after  substi- 
,  AL 

tution  of  —  , 
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A3 


(4.6) 


W2a(«) 


+  Wyc($)  + 


2 

w  e($) 


-1 


where 


W 


(4.7) 


n 


(4.8) 


u 


(4.9) 


and  a(4>),  c(<j>),  and  e(<}>)  are  functions  of  A. 


a(<(>) 


27 
4  <!>  ^ 


Jl  _  + _ 1  |  3cot4> 

6  "  2  .  ,  2^  +  4$ 

<(>  4sin  4> 


(4.10) 


c(0)  H  -2L- 

40>4 


JL _  _  cot4> 


U2  2sin2<*>  21,1 


(4.11) 


e  (<(>)  = 


4<t2  L 


I  +  _L__ 
6  ,  4  2 
4sm  ij) 


cotft 

4<j> 


(4.12) 


In  arriving  at  this  result,  power  series  in  (M  have  been  truncated; 
for  example, 

a(<j>)  =  a^*)  +  (|)  a2(t)  + . 


has  been  truncated,  retaining  only  the  first  term,  which  is  consistent 
with  the  assumption  of  small  curvature. 

As  $-K),  the  above  expressions  lose  numerical  accuracy  badly.  For 
small  values  of  <(>,  the  following  expressions  should  be  used. 


44 


45 


FIGURE  4  2  DEFLECTION  BEHAVIOR  OF  CURVED  BEA 


Four  curves  are  shown  in  Figure  4.2,  corresponding  to  four  differ¬ 
ent  values  of  u.  For  any  value  of  p,  <f=0  gives  ui=0  and  w=7;  hence,  the 
ends  of  the  curves  fall  at  these  two  values  of  u>.  Consider  first  the 
curve  for  p=l,  a  very  small  value  of  p  corresponding  to  a  nearly  straight 
beam.  As  cp  is  increased  from  zero,  the  two  values  of  w  found  for  each 
value  of  <j>  move  from  the  initial  values,  0  and  7,  toward  each  other,  and 
meet  at  a  small  value  of  <p .  Higher  values  of  <J>  would  give  a  complex 
solution  for  w.  For  small  values  of  p,  the  normalized  deflection  is 
essentially  linear  with  loading. 

Next  consider  the  deflection  behavior  for  a  moderately  large  value 
of  it ,  p  -^/lOOO.  As  $  increases  from  zero,  the  two  values  of  uj  approach 
one  another,  and  are  equal  at  p=n.  The  value  of  ui  corresponding  to  c=Tt 
is  o)=7t2/3=3. 290.  As  <p  increases  beyond  it,  the  two  branches  behave  as 
shown,  eventually  converging  at  a  point  slightly  above  the  lower  knee  of 
a  value  of  tp  of  a  little  over  5.  Having  established  the  curve,  consider 
the  structural  behavior  as  the  loading  parameter,  u,  increases.  Start¬ 
ing  at  oi=0,  the  axial  force  parameter,  $>,  is  zero.  As  oj  increases,  the 
deflection  increases  and  so  does  4> .  When  to  reaches  the  critical  value 
defined  by  the  lower  knee,  a  further  increase  in  w  must  be  accompanied 
by  a  jump  in  deflection,  as  shown  by  the  dashed  line.  At  the  same  time, 
there  is  a  step  decrease  in  $.  This  is  the  snap-through  phenomenon.  At 
the  lower  knee,  the  deflection  pattern  leaves  the  center  of  the  beam 
above  the  chord  joining  the  two  ends  of  the  beam;  at  the  end  of  snap- 
through,  the  center  of  the  beam  is  well  below  this  chord. 

Curves  for  two  other  values  of  p  are  also  shown  in  Figure  4.2.  The 
first,  for  p=°°,  exhibits  the  same  general  behavior  as  the  curve  for 
p  -^1000.  In  this  case,  however,  there  is  no  deflection  until  the 
snap-through  load,  4tt 2 / 3 ,  is  reached. 

The  final  curve  is  for  p  =  2tt2 which  is  the  critical  value  of 

p,  pcr.  For  p<pcr,  snap-through  does  not  occur  (that  is,  the  curve  is 

single-valued  for  all  values  of  w);  for  p>p  ,  snap-through  does  occur. 

For  p=p  ,  the  maximum  value  of  6  is  P=tt. 
cr 

Next  consider  the  behavior  of  the  beam  beyond  the  terminal  point  of 
the  upper  branch  on  Figure  4.2.  As  the  terminal  point  is  approached,  is 
decreasing  and  reaches  zero  at  the  terminal  point.  This  means  that  the 
compressive  axial  load  decreases  to  zero  at  the  terminal  point.  For 
higher  loadings  and  deflections,  the  axial  load  becomes  negative;  that 
is,  the  axial  load  becomes  tensile.  The  loading  parameter,  $,  then 
becomes  imaginary,  as  may  be  seen  from  Eq.  (4.4).  For  convenience,  the 
loading  parameter  £  is  used  in  the  tensile  regime,  where 


<f>  =  U 


(4.19) 


Referring  again  to  Figure  4.2,  the  following  types  of  elastic  behavior 
are  possible: 
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1) 


y  <  u 


cr 

Snap-through  does  not  occur 


2) 


V  >  V 


cr 


a) 


w  at  snap-through  <7 

Beam  snaps  through  into  compression  regime 
w  at  snap-through  >7 

Beam  snaps  through  into  tension  regime 


In  order  to  interpret  beam  response  in  terms  of  damage,  stresses 
must  be  found.  Since  the  deformation  shape  has  been  determined,  this 
may  easily  be  accomplished,  resulting  in 


n2  |  =  -  f-  +  f>\  (l-<f>cot<t>) 
for  a  compressive  axial  force  and 


(;M 


(4.20) 


2  a 

n2  = 


E  3  2 

for  a  tensile  •  ial  force,  where 


(CcothC-1) 


(?-) 


(4.21) 


2c 

x  =  rn 


(4.22) 


Here  c  is  the  distance  from  the  beam  neutral  axis  to  the  outer  fiber  on 
the  loaded  side  of  the  beam.  The  stresses  given  by  Eqs.  (4.20)  and 
(4.21)  are  those  which  occur  at  the  clamped  support  and  are  the  maximum 
stresses  which  occur.  In  addition,  Eqs.  (4.20)  and  (4.21)  are  chosen  to 
provide  the  tensile  stresses  at  the  support,  rather  than  the  compressive 
stresses,  since  the  vulnerability  program  will  be  directed  toward  ten¬ 
sile  failure.  In  passing,  it  is  noted  that  the  stress  which  occurs  when 
the  axial  stress  switches  from  compression  to  tension  is 


=  3.5pX  (4.23) 

<j>=0,  upper  branch 
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The  parameters  n  and  x  which  have  been  defined  in  the  present 
subsection  are  functions  of  the  beam  cross-section.  Some  insight  into 
what  the  parameters  represent  is  attained  by  noting  that  for  a  solid 
cross-section  n  is  the  beam  length  to  thickness  ratio  and  x  is  unity. 

4.3  Straight  Beam  Results 

Reduction  of  the  curved-beam  solution  of  the  foregoing  subsection 
to  the  straight-beam  case  is  fairly  straightforward.  Substituting  p=0 
into  Eq.  (4.6)  and  recognizing  that  the  axial  force  will  always  be 
tensile  for  a  straight  beam, 


W  =  iJ-l/aU) 


(4.24) 


where 


a(0  = 


27 

4£f 


2 

4sinh  £ 


3coth£ 

4£ 


t>0.2 


(4.25) 


a(£)  =  - 


140£ 


+  0.00143,  £<0.2 


(4.26) 


The  stress  relation  for  the  straight  beam  may  be  derived  from  Eq.  (4.21) 
as 


n*  I  -  f  +  !*•  isasfci  „  .  (4.27) 

4.4  Fitting  Functions  for  2024  A1  Straight-Beam  Data 

As  indicated  in  Subsection  4.1,  the  straight-beam  data  were  fitted 
first,  and  adjustments  were  made  to  the  straight-beam  fits  to  accommo¬ 
date  the  curved  beams.  Futhermore,  for  the  straight  beams,  the  light 
damage  data  were  fitted  first,  since  it  seemed  likely  that  guidance 
provided  by  the  static  elastic  solution  would  be  most  appropriate  for 
the  light-damage  case.  Hence,  the  straight-beam  fitting  functions  for 
light  damage  are  presented  first. 

4.4.1  Light  Damage  Fitting  Functions 

The  relevant  static  result  is  Eq.  (4.27),  which  must  be  solved 
together  with  Eqs.  (4.24)  and  (4.25)  or  (4.26).  Equation  (4.27)  yields 
the  elastic  strain  resulting  from  the  application  of  a  static,  uniform 
pressure.  If  a  step-function  pressure  (corresponding  to  the  pressure 
asymptote  on  the  P-I  hyperbola)  were  applied  to  the  structure,  a  dynamic 
overstress  factor  would  be  expected  compared  with  the  stress  or  strain 
obtained  by  static  application  of  the  same  pressure.  Hence,  one  might 
expect  that  Eq.  (4.27)  would  apply  approximately  for  a  step-function 
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pressure  if  a  dynamic  factor  (really,  the  inverse  of  the  dynamic  factor) 
were  applied  to  the  left-hand  side  of  Eq.  (4.27).  This  was  attempted, 
and  it  was  found  that  excellent  correlation  with  the  computer  results 
could  be  obtained  by  making  this  factor  a  function  of  x-  The  modified 
version  of  Eq.  (4.27)  used  to  solve  for  the  pressure  asymptote  thus 
becomes 


n2s(x)  =  J-  +  ^  ;c°t2~1  w  (4.28) 

4 

Solution  of  this  equation  in  conjunction  with  Eqs.  (4.24)  and  (4.25)  or 
(4.26)  gives  the  pressure  asymptote,  P  .  The  function  s(x)  is  plotted 
in  Fig.  4.3.  0 

Having  successfully  fitted  the  P0  data,  attention  is  next 
focused  on  the  impulse  asymptote,  I0.  The  approach  used  is  to  relate 
the  impulse  to  the  work  done  by  the  application  of  P  ,  which  is  given 
by  ° 


(4.29) 


where  f(4)  is  a 


'inction  of  4  which  results  from  the  indicated  integration. 


The  kinetic  energy  imparted  by  the  impulse  is 


KE  =  y  pALV2 


(4.30) 


where  V  is  the  velocity  resulting  from  the  impulse  and  is  given  by 


V 


(4.31) 


Substituting  Eq.  (4.31)  into  Eq.  (4.30)  and  equating  the  kinetic  energy 
imparted  by  the  impulse  to  the  work  done  by  Po, 


2P  2L4 
o 


El 


f(4) 


(4.32) 
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FIGURE 


Introducing  '  as  defined  by  Eq .  (4.8),  and  defining  W  consistently  with 
Eq.  (4.7)  as  ° 


P  I. 

O  ; 

AE 


(4.33) 


Eq .  (4 .  32  )  be i  ones 


1 

(4.34) 

y.y\ 

where  '(  )  is  simplv  a  constant  times  the  square  root  of  f(0*  Equation 

(4.34)  may  be  used  to  find  l  ;  however,  it  was  found  that  better  corre¬ 
lation  with  the  computer  data  was  obtained  by  introducing  a  function  of 
the  parameter  2c/t. 


where 


a 


1.71 


(4.35) 


(4.36) 


and  <5(0  is  plotted  in  Figure  4.4.  The  parameter  c  is  the  value  of  4 
corresponding  to  the  solution  of  Eq.  (4.28).  It  should  be  noted  that 
<5(0  was  actually  determined  empirically  from  the  computer  data  rather 
than  through  integration  of  the  deformation  shape. 


The  final  step  is  to  establish  a  fitting  function  for  the 
constant  J  in  Eq.  (3.22).  It  may  be  noted  in  passing  that  accurate  fitting 
of  J  is  not  nearly  so  critical  as  the  fitting  of  PQ  and  I0.  For  example, 

PI 

at  the  knee  (defined  by  *  y) 
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r  -  1 

o 

The  constant  J  falls  typically  in  the  range  of  0.4  to  1.0.  As  an 
example,  with  J  at  an  average  value  of  0.7,  a  fitted  value  of  1.0  for  .1, 
which  represents  over  a  40%  error  in  fitting,  would  result  in  only  a  9". 
error  ia  P  . 

Fitting  of  the  constant  J  was  done  purely  on  the  basis  of  the 
computer  data,  resulting  in  the  fitting  function 


J  =  j(x)6(0 


(4.37) 


where 


j(X)  i  0.57  +  0.52X 


(4.38) 


and  6(4)  is  plotted  in  Figure  4.5. 

4.4.2  Moderate  and  Heavy  Damage  Fitting  Functions 

The  pfocedure  used  to  fit  the  computer  data  for  moderate  and 
heavy  damage  was  to  build  on  the  fitting  functions  for  light  damage.  v 
The  factors  used  to  convert  the  light  damage  factors  to  moderate  and 
heavy  damage  factors  were  derived  empirically  from  the  computer  data. 

It  was  found  that  factors  could  be  based  upon  4,  which,  as  indicated 
earlier,  is  the  value  of  4  corresponding  to  the  solution  of  Eq.  (4.28) 

P  P 

°2  °3 

for  light  damage.  Figure  4.6  defines  the  ratios  - —  and  - — ,  which 

°1  °1 

are  the  ratios  of  P0  for  moderate  and  heavy  damage,  respectively,  to  the 
value  of  PQ  for  light  damage.  Figure  4.7  presents  similar  results  for 
the  impulse  asymptote,  Iq. 

For  the  constant  J,  Eqs.  (4.37)  and  (4.38)  are  used,  but 
8(4)  is  redefined  in  Figure  4.8  for  moderate  and  heavy  damage. 

4.5  Fitting  Functions  for  2024  A1  Curved-Beam  Data 

As  pointed  out  earlier,  the  presence  of  the  snap-through  phenomenon 
greatly  complicates  the  curved-beam  solution  in  comparison  with  the 
straight -beam  solution.  Using  the  analysis  in  Subsection  4.1,  the 
static  value  of  «  which  produces  snap-through  can  be  calculated.  This 

g 

value  of  u  will  be  designated  by  w  ,  where  the  subscript  indicates 

O  1 

snap-through  and  the  superscript  denotes  static.  Also,  using  Eqs. 
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<  5)6/ 
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FIGURE  4  5/3(0  VS  l  FOR  LIGHT  DAMAGE 


(A. 20)  and  (A. 21),  the  values  of  rp  —  just  before  and  just  after  snap- 

through  can  be  established.  These  will  be  designated  by  SN2B  and  SN2A, 

respectively  where  the  notation  SN2  is  used  to  repalce  n2  ^  and  the  B 
and  A  indicate  before  and  after  snap-through. 

In  some  cases,  stress  or  strain,  rather  than  snap-through,  will  be 
the  dominant  factor.  From  the  straight-beam  fitting  functions,  the 

parameter  s(y)  is  available  which  defines  an  effective  value  of  ~. 

E 

The  quantity  n2s(x)  will  be  denoted  by  SN2.  The  static  value  of  u>  will 

g 

be  denoted  by  w  ,  where  the  subscript  indicates  that  uj  is  based  upon 
stress  considerations. 

In  many  cases,  SN2  will  occur  in  the  snap-through  region;  that  is, 
SN2B£SN2<SN2A.  In  such  cases,  will  be  set  equal  to  w®  . 

s 

If  U<Ucr,  there  is  no  snap-through  and  hence  will  be  undefined. 

In  the  procedure  to  be  used  for  fitting  the  curved  beam  data,  this 
situation  would  lead  to  discontinuities  in  the  solution.  To  avoid  such 
discontinuities,  an  effective  snap-through  loading,  u>,  is  defined  as 
follows.  The  maximum  value  of  the  axial  load  parameter,  <J> ,  is  found. 

From  Subsection  A. 2,  it  is  known  that  4>max<’T  for  p<pcr.  Corresponding 

S  9  O 

to  <p  are  a  value  of  w,  which  is  taken  as  cu  ,  and  a  value  of  rp  — ,  which 
max  SI  L 

is  used  to  define  both  SN2B  and  SN2A. 

A. 5.1  Light  Damage  Fitting  Functions 

The  various  parameters  developed  in  the  previous  section  are 
used  to  define  the  light  damage  fitting  functions.  In  fitting  PD, 

s  s  s  s 

there  are  two  main  cases;  u>  <  <*>_  and  u  >  u)„_.  These  two  cases  are 
,,,  o  Sjl  o  —  ST 

examined  below. 

s  s 

Case  1:  u  < 

_ 0 _ ST 

The  desired  value  of  u>,  from  which  P0  can  be  calculated, 
is  the  smaller  of  or  where  the  superscript  denotes  dynamic  and 


|  1.0, 

I  1.0-0.A5  -  O.aJ, 


|<  0.A 


>  0. A 


(4.39) 


(pcr’  k) 

("■  i) 


"  i  “cr 


P  >  u 


(A. AO) 
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In  Eq.  (A. 40),  Fc 


(“*  R)l 


is  a  snap-through  factor  which  is  plotted  in 


Fig.  4.9.  The  curve  for  j  =  0  is  shown  only  to  permit  interpolation  of 
L  K 

F  for  values  of  —  less  than  0.4;  it  is  apparent  that  p>0  is  not  com- 
b  i  ^  R 

patible  with  —  =  0. 

In  preparation  for  the  fitting  of  the  moderate  and  heavy  damage 
computer  data,  SN2B  is  modified  as  follows: 


SN2B  = 


uiD  -  CJD 

SN2  +  £  (SN2B-SN2) ,  wjL  >  uP 

s  D  ST  o 

U5T  -  U> 

ST  a 


D  D 

UCT  * 

ST  —  o 


(4.41) 


s  s 

Case  2:  oj  < 

_ o  —  ST 

For  this  case,  the  desired  value  of  m  is  approximated  by 


«D  ,  I  D  D  \  J  Ucr 
"ST  +  (“o  "  “st)  V  ~r  * 


p  >  p 


(4.42) 


where 


,  u  1  u 


M) 

(“«•  s) 


»  V  >  V, 


,  M  £  Vr 


(4.43) 


M) 

('«•  s) 


,  V  >  W, 


,  n  in. 


(4.44) 
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FIGURE  4.9  SNAP-THROUGH  FACTOR,  FgT 


As  in  the  case  of  the  straight  beam,  the  curved-beam  fitting 
function  for  the  impulse  asymptote,  IQ,  is  based  upon  energy  consider¬ 
ations-  Basically,  Eq.  (4.35)  is  used  for  the  curved  beam,  with  two 
modifications.  The  first  has  to  do  with  the  deflection  function,  6.  It 
will  be  recalled  that  the  value  of  6  used  for  the  straight  beam  was  that 
determined  empirically  from  the  computer  data,  rather  than  the  value 
calculated  by  integration  of  the  deflection  shape.  It  is  apparent  that 
the  straight-beam  5  is  inappropriate  for  curved  beams.  The  procedure 
adapted  was  to  determine  6  analytically,  from  the  deflection_shape ,  for 
both  the  straight  and  curved  beams.  These  are  denoted  by  6(0  and  6(0, 
respectively,  where  <j>  is  the  value  of  <f>  corresponding  to  the  solution  of 
Eq.  (4.6)  and  4  is,  as  indicated  earlier,  the  value  of  C  corresponding 
to  the  solution  of  Eq.  (4.28)  for  the  straight  beam.  In  the  case  in 
which  the  curved-beam  solution  lies  in  the  snap-through  region,  6(0  cannot 
be  defined  directly;  however,  the  values  of  6  just  before  and  just  after 
snap-through  can  be  found.  These  will  be  denoted  by  6^  and  6^,  respec¬ 
tively.  The  equivalent  value  of  6  in  the  snap-through  region  is  cal¬ 
culated  from 


6(0 


SN2-SN2B  -2\ 

SN2A-SN2B  (6A  6B  ) 


(4.45) 


The  expression  which  defines  6(0  based  on  integration  of  the 
deflection  shape  is 


6(0 


(4.46) 


Equation  (4.46)  is  used  by  itself  to  define  6  outside  of  the  snap- 
throughregion^  In  the  snap-through  region,  Eq.  (4.46)  is  used  to 
define  6fi  and  6^,  and  Eq.  (4.45)  is  then  used  to  determine  6(0. 

The  expression  for  6(c)  for  the  straight  beam  is 


6(C) 


-  ccothc  + 


?)] 


1/2 


(4.47) 


This  equation  is  poorly  behaved  for  small  values  of  C»  as  is  Eq.  (4.46) 
for  small  values  of  <j>.  The  relevant  expressions  for  small  values  of 
$  and  c  are 
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$  1  0.3 


(4.48) 


4  <  0.3  (4.49) 


The  modification  to  Eq.  (4.35)  to  reflect  the  curved-beam 


deflection  is  simply  to  replace  6(£)  by  6(<p) 


It  is  clear  that 


_6(0 

this  reduces  to  the  straight-beam  result  for  very  small  curvature,  since 
for  that  case  6($)  'v.  &(r).  In  this  connection,  it  must  be  noted  that  $ 
for  this  case  will  be  imaginary,  indicating  axial  tension.  The  use  of 
4>  instead  of  i£  or  simply  c  is  retained  to  simplify  differentiation 
between  the  curved  and  straight  beam. 


The 

cation  by  the 


second  modification  made  to  Eq.  (4.35)  involves  multipli- 
factor  used  in  Eq.  (4.39).  Thus,  the  final  result  is 


(4 


The  fitting  function  for  the  constant  J  involves  a  very  simple 
modification  of  Eq.  (4.37). 


J  -  1(X)6(0  (l  + 


(4.51) 


4.5.2  Moderate  and  Heavy  Damage  Fitting  Functions 

P 

o 

From  the  straight-beam  solution,  the  parameter  ^ -  is  avail¬ 
able.  A  curved-beam  pressure  factor  is  defined  as  o 


(4.52) 


50) 
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From  the  curved-beam  solution  for  light  damage,  the  loading  parameters 

U|J  ,  and  are  available,  where  the  subscript  "1"  has  been  added 

to  denote  the  first  damage  level,  light  damage.  A  new  loading  parameter 
is  defined  by 


1 

( 


SN2A 

Max(SN2, SN2B) 


UJ 


D 

ST 


SN2A>SN2 


SN2A<=SN2 


(4.53) 


where  the  notation  Max(a,b)  signifies  that  the  larger  of  "a"  or  b 

should  be  used.  The  loading  range  from  to  may  be  considered 

to  represent  an  effective  dynamic  loading  range  for  snap-through. 
Finally,  one  more  loading  parameter  is  defined 


u>  H  u^PF  (4.54) 

n  1  n 

Based  on  the  definitions  given  by  Eqs.  (4.52)  -  (4.54),  the  desired 
loading  parameter  which  gives  the  pressure  asymptote  for  moderate  and 
heavy  damage  may  be  found  as  follows: 

/ 

Case  1:  w  <  u>crr 
n  —  ST 


w  *  u>  (4.55) 

n  n 


Case  2: 


0} 

n 


(4.56) 


Case  3: 


w  > 

n 


Define  another  Intermediate  loading  parameter  by 


(0 

n 


w  -wc> 
n  of 


D 

ST 


(4.57) 
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Then 


(4.58) 


Turning  next  to 
straight-beam  values 


the  impulse 
P 

o 


of  -  and 


asymptote,  use  will  be  made  of  the 
I 

o 

_ n 

I 


o 


1 


(4.59) 


where  the  subscript  "CB"  denotes  curved-beam  and  is  used  to  differen¬ 
tiate  the  curved-beam  impulse  ratio  from  the  straight-beam  impulse  ratio 
on  the  right-hand  side  of  Eq.  (4.59). 

Finally,  the  constant  J  for  the  curved  beam  for  moderate  and 
heavy  damage  is  given  by  Eq.  (4.51),  with  the  value  of  8(0  selected  for 
the  appropriate  damage  level. 

4.6  Application  of  Fitting  Functions  to  7075  A1  and  to  Plates 

The  fitting  functions  defined  thus  far  apply  to  straight  and  curved 
beams  of  2024  Al.  For  7075  A1  beams,  the  pressure  and  impulse  asymptote 
may  be  found  by  multiplying  the  2024  Al  results  by  the  following  fac¬ 
tors;  1.4,  0.95,  and  0.86  for  the  pressure  asymptote,  and  1.3,  0.71,  and 
0.68  for  the  impulse  asymptote,  where  the  three  numbers  apply  for  the 
three  damage  levels.  In  order  to  avoid  situations  in  which,  for  example, 
P  <  P  ,  which  might  occur  in  snap-through  cases  using  the  above 
°2  °1 

,  numbers,  P  is  taken  as  the  larger  of  P  and  P  ,  etc.  The  constant  J 
o„  o,  o_ 

\  2  12 

N  is  unchanged. 

For  plates,  the  beam  results  are  multiplied  by  a  function  of  aspect 
ratio,  AR,  where  the  aspect  ratio  is  the  ratio  of  the  length  of  the  long 
side  of  the  plate  to  the  length  of  the  short  side.  The  beam  is  effec¬ 
tively  a  plate  of  infinite  aspect  ratio.  The  function  of  aspect  ratio 
is 
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-ll(AR-l) 


(4.60) 


F 


AR 


1  +  .47e 


Both  the  pressure  and  impulse  asymptotes  for  the  beam  are  to  be  mul¬ 
tiplied  by  F  .  The  constant  J  is  unchanged. 

4.7  Fitting  Accuracy 

The  estimated  accuracy  of  the  fits  defined  in  the  preceding  sub¬ 
sections  is  given  in  Table  4.1.  The  accuracies  cited  should  be  regarded 
as  seldom-exceeded  values;  that  is,  in  most  cases  the  accuracy  will  be 
better  than  that  given.  The  fitting  functions  are  believed  to  be 
unbiased;  that  is,  on  the  average  the  fitting  functions  will  agree  with 
the  computer  data.  It  is  important  to  realize  that  the  accuracies  pre¬ 
sented  describe  the  fitting  accuracy  only;  they  do  not  purport  to  define 
the  overall  accuracy  of  the  procedure  in  predicting  various  levels  of 
damage  for  a  real  structure.  The  latter  could  be  realistically  defined 
only  through  a  comprehensive  testing  program  and/or  comparison  with 
existing  test  data. 
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TABLE  4.1 

ACCURACY  OF  FITTING  FUNCTIONS 


Structural  Type 

Straight 

Beam 

Curved 

Beam 

Flat 

Plate 

Curved 

Plate 

Material 

2024 

7075 

2024 

7075 

2024 

7075 

2024 

7075 

Estimated 
Accuracy  (%) 

5 

10 

15 

20 

15 

20 

25 

25 
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SECTION  5 

PIVUL  PROGRAM  DESCRIPTION 


5 . 1  Introduction 

This  section  describes  PIVUL,  the  Pressure-Impulse  VULnerability 
program.  The  program  utilizes  the  data  from  the  computer  runs,  fitted 
as  described  in  the  previous  section,  to  predict  vulnerability  of  simple 
structural  elements.  There  are  four  basic  vulnerability  options: 

1.  The  user  specifies  the  weapon  yield,  defines  the  struc¬ 
ture  and  the  desired  level  of  damage,  and  the  program 
calculates  the  range  at  which  that  level  of  damage 
occurs . 

2.  The  user  specifies  the  duration  (tQ)  of  an  exponential 

t 

forcing  function  (P  =  PM^  |l  -  )e  to  )  ,  defines  the 

o 

structure  and  the  desired  level  of  damage,  and  the  program 
calculates  the  value  of  Pw>„  which  produces  that  level 
of  damage. 

3.  The  user  specifies  the  weapon  yield  and  range  and  defines 
the  structure  and  the  program  calculates  the  correspond¬ 
ing  level  of  damage. 

4.  The  user  provides  a  curve  of  pressure  vs  time  and  defines 
the  structure  and  the  program  calculates  the  level  of 
damage . 

In  all  cases,  ambient  conditions  are  assumed  to  correspond  to  sea  level. 
Options  1  and  2  are  solved  by  trial  and  error  by  the  program;  the  solu¬ 
tions  for  Options  3  and  4  are  direct. 

As  indicated  earlier,  the  loading  on  the  cylinder  is  assumed  to  be 
sufficiently  long  in  comparison  with  the  response  time  so  that  only  the 
peak  overpressure  is  required  to  define  the  loading.  Hence,  Options  2 
and  4  are  inapplicable  for  cylinders. 

For  plates  and  beams,  the  surface  of  the  structure  containing  the 
plate  or  beam  is  assumed,  for  Options  1-3,  to  be  oriented  parallel  to 
the  shock  front,  since  this  is  the  most  critical  orientation,  or  at 
least  very  close  to  it.  The  loading  produced  by  the  shock  reflection 
and  subsequent  diffraction  waves  is  included  in  the  loading  calculated 
by  the  program  for  Options  1-3. 

Subsection  5.2  contains  a  brief  description  of  the  various  sub¬ 
routines  contained  in  PIVUL.  The  major  program  variables  are  defined 
in  Subsection  5.3  and  the  program  input  is  described  in  Subsection  5.4. 
Subsection  5.5  contains  comments  on  program  operation  and  describes  the 
output.  Finally,  an  example  problem  is  given  in  Subsection  5.6. 
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PRECEDING  PAOB  BLANK-NOT  FILMED 


5.2  Description  of  Subroutines 

The  PIVUL  program  is  the  main  program  which  controls  the  flow  of 
the  overall  program.  PIVUL  also  reads  the  input  and  provides  the 
output.  PIVUL  calls  the  following  subroutines:  BEAM,  BLAST,  CYL, 

PCURV ,  PHYP,  PINT,  RITER,  SIDE,  YOUNG  and  YPI. 

Of  the  above  group  of  subroutines,  only  BEAM,  BLAST,  and  CYL  call 
other  subroutines.  The  order  of  listing  the  subroutines  will  fall  into 
three  groups  below;  those  called  by  PIVUL,  except  for  BEAM,  BLAST,  and 
CYL,  the  BEAM  group,  the  BLAST  group,  and  the  CYL  group.  Within  each 
group,  the  order  will  be  alphabetical  following  the  master  program. 

5.2.1  PIVUL  Group 


PCURV 

Given  peak  pressure  and  duration  of  an  exponential  forcing  function, 
calculates  pressure-time  points. 

Called  by  PIVUL. 


PHYP 

Calculates  intersection  of  straight  line  with  hyperbola  in  P-I  plane. 
Called  by  PIVUL. 


PINT 

Calculates  pressure  and  impulse  from  points  defining  pressure-time 
wave . 

Called  by  PIVUL. 


RITER 

Controls  iteration  to  determine  critical  range  for  option  1  or  critical 
pressure  for  option  2. 

Called  by  PIVUL. 


SIDE 

Calculates  side-on  loading  from  free-field  pressure-time  history. 
Called  by  PIVUL. 
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YOUNG 


Calculates  effective  pressure  and  effective  impulse  from  points  defining 
pressure-time  curve  using  modified  Youngdahl  approach. 

Called  by  PIVUL. 


YPI 

Defines  the  exponential  forcing  function  which  gives  desired  effective 
pressure  and  effective  impulse  and  calculates  associated  pressure  and 
impulse . 

Called  by  PIVUL. 

5.2.2  BEAM  Group 


BEAM 

Master  program  for  BEAM  group.  Defines  P-I  hyperbola  for  straight  beams 
and  controls  flow  for  curved  beams. 

Calls  BDATA,  CURVB,  INTI,  and  SNAP. 

Called  by  PIVUL. 


BDATA 

Contains  data  defining  fitting  functions  described  in  Section  4. 
Called  by  BEAM. 


CURVB 

Defines  P-I  hyperbola  for  curved  beams. 
Calls  INTI,  INTERP  and  SETUP. 

Called  by  BEAM. 


INTI 

Performs  one-dimensional  interpolation. 
Called  by  BEAM  and  CURVE. 
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INTERP 


Performs  one-,  two-,  and  three-dimensional  interpolations. 
Called  by  CURVB  (and  by  CYL  in  CYL  group). 

SETUP 

Sets  up  tables  of  functions  of  <p  used  in  solving  Eq.  (4.17). 
Called  by  CURVB. 


SNAP 

Calculates  value  of  oj  associated  with  snap-through  and  value  of  u; 
associated  with  the  stress  parameter  SN2.  Provides  values  of  stress 
parameter  before  and  after  snap-through,  and  associated  deflection 
parameters  for  use  in  defining  the  impulse  asymptote  on  the  P-I 
hyperbola. 

Calls  WTEN. 

Called  by  REAM. 


WTEN 

Calculates  w  for  solution  in  tension  region. 
Called  by  SNAP. 

5.2.3  BLAST  Group 


BLAST 

Calculates  pressure-time  points  given  weapon  yield  and  range. 
Calls  HYPRA  and  WFZR. 

Called  by  PIVUL. 


HYPRA 

Processes  AFWL  1KT  standard  blast  model. 
Calls  I0PT1,  I0PT2 ,  and  I0PT3. 

Called  by  BLAST. 
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I0PT1 


Determines  overpressure  given  time  and  range. 
Calls  KSHK. 

Called  by  HYPRA. 


I0PT2 

\ 

Determines  blast  wave  radius  and  peak  overpressure  given  time  (not  used  \ 
in  PIVUL  program) . 

Calls  WFPKOP  and  WFPR.  \ 

Called  by  HYPRA. 

I0PT3 

Determines  time  of  shock  arrival  and  peak  overpressure  given  range. 

Calls  WFPKOP  and  WFPR. 

Called  by  HYPRA. 

KSHK* 

Calculates  overpressure  given  time  and  range. 

Calls  WFPKOP,  WFPR,  WFPRMT ,  and  WFZR. 

Called  by  I0PT1. 


WFPKOP* 

Calculates  peak  overpressure  at  specified  range. 
Called  by  I0PT2,  I0PT3,  and  KSHK. 


WFPR* 

Calculates  blast  wave  radius  at  specified  time. 
Called  by  I0PT2  and  KSHK. 


These  routines  are  modifications  of  routines  written  by  AFWL. 
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WFPRMT* 


Calculates  overpressure  waveform. 
Called  by  KSHK. 


WFZR* 

Calculates  the  range  corresponding  to  zero  overpressure  at  specified 
time. 

Called  by  BLAST  and  KSHK. 

5.2.4  CYL  Group 


CYL 

Determines  the  critical  free-field  overpressure  for  a  cylinder  from 
tabulated  data.  The  tabulated  data  are  functions  of  R/t,  L/D,  and 
damage  level. 

Calls  INTERP. 

Called  by  PIVUL. 


INTERP 

Performs  one-,  two-,  and  three-dimensional  interpolations. 
Called  by  CYL  (and  by  CURVB  in  BEAM  group). 


5 . 3  Major  Program  Variables 

The  major  program  variables  are  defined  in  the  following  list. 

An  asterisk  preceding  a  variable  name  indicates  that  the  variable  is 
an  input  variable.  The  dimension  of  the  variable  is  also  given, 
if  applicable.  Most  variables  are  in  English  units,  since  the  program 
works  in  that  system.  All  input-output  quantities  are  metric.  Varia¬ 
bles  which  are  metric  end  in  the  letter  M,  and  the  following  list 
specifically  identifies  metric  variables. 
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A(629) 


a ( (p )  for  629  values  of  $>,  where  a(p)  is  defined  in 
Eqs .  (4.10)  and  (4.13) 


AEL  Cross-sectional  area  of  beam  times  modulus  of  el^^ixcity 

divided  by  beam  length 

AJ(3)  Hyperbolic  constants,  J,  for  three  damage  levels 

ALPHA  a  defined  in  Eq.  (4.36) 


AR 

ATOTL 

*AT0TM 

BETA 

BETAT (186) 

BI 

BIGW 

*BLM 

BLNTH 

BP 

CAE(629) 

CENTI 

CHI 

CHIT(31) 

*CLDD 
*CLDR 
CLDRT (3) 

CNST 


Plate  aspect  ratio 

Cross-section  area 

Cross-section  area,  metric 

Value  of  8(4)  used  to  define  J,  Eq .  (4.37) 

Tabulated  values  of  fitting  function  6(4)  plotted  in 
Figures  4.5  and  4.8 

Impulse  or  effective  impulse  from  blast  pulse 
Loading  parameter,  W,  defined  in  Eq.  (4.7) 

Beam  length,  metric 
Beam  length 

Pressure  or  effective  pressure  from  blast  pulse 

for  629  values  of  <p ,  where  a  (<}>),  c($) 
Eqs.  (4.10  -  4.15) 

Cross-section  moment  of  inertia 

X,  defined  in  Eq.  (4.22) 

Values  of  x  used  in  tabulating  data  for  s(x)  in  Figure 
4.3 

Cylinder  length  to  diameter  ratio 

Beam  length  to  radius  of  curvature  ratio 

Values  of  L/R  used  in  tabulating  data  for  snap-through 
factor  plotted  in  Figure  4.9 

Hyperbolic  constant  J 


C2(t) 

-  ew>J 

are  defined  in 
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*CNi'IM 

CNTRD 

*CNTRM 

CRIT (4) 

C2A(629) 

DBA 

DBAR 

DRB 

DELT(31) 

DELZ 
EMOD 
ET  (2) 

ETA 

FACT 

FACTI (90) 

FIM(3) 

FPM(3) 


Cross-section  area  moment  of  inertia,  metric 

Distance  from  neutral  axis  to  outer  fiber  on  loaded 
side 

Distance  from  neutral  axis  to  outer  fiber  on  loaded 
side,  metric 

Response  criteria  used  in  iteration  process  which  drives 
criterion  to  1.0  in  CRIT (N) ,N=1 , 2 , 3;  CRIT(4)  contains 
convergence  criterion  for  iteration 

for  629  values  of  $,  where  a (4> )  and  c(i}>)  are 

2a(<J>) 

defined  in  Eqs.  (4.10),  (4.11),  (4.13),  and  (4.14) 
Deflection  parameter  $A  which  appears  in  Eq.  (4.45) 
Deflection  parameter  6  defined  in  Eqs.  (4.45)  -  (4.49)) 
Deflection  parameter  6  which  appears  in  Eq.  (4.45) 

O 

Tabulated  values  of  fitting  function  6(5)  plotted  in 
Figure  4.4 

Value  of  5(5)  used  to  define  Iq,  Eq.  (4.35) 

Modulus  of  elasticity 

Table  of  moduli  of  elasticity 

n,  defined  in  Eq.  (4.8) 

Value  of  I  /I 

o  o. 
n  1 

Tabulated  values  of  fitting  function  I  /I  ,  plotted 
in  Figure  4.7  n  1 

Impulse  factor  for  7075-T6A1  compared  with  2024-T4A1 

Pressure  asymptote  factor  for  7075-T6A1  compared  with 
2024-T4A1 


FST 


Snap-through  factor 


*IBPC 


Structure  code 

1- beam 

2- panel 

3- cylinder 
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I 


I 

i 

I  *IMAT 


ISC (90) 
*ISL 


*IVOPT 


KOK 

NC 

NCHI 

NLDR 

*NPP 

NTRIAL 

NU 

NZ 

NZETA 
NZTAI 
NZTA1 
OUT (2) 

PCR 

PCRM 

PFACT 


Material  code 

1- 2024T4  Aluminum 

2- 7075T6  Aluminum 

Output  array  used  to  indicate  level  of  damage 

Damage  level  code 

1- light  damage 

2- moderate  damage 

3- heavy  damage 

Vulnerability  option  code 

1- given  weapon  yield,  determine  critical  range 

2- given  duration  of  exponential  forcing  function, 
determine  critical  pressure  level 

3- given  weapon  yield  and  range,  determine  damage 
level 

4- given  curve  of  pressure  vs  time,  determine  damage 
level 

Control  constant  set  to  1  by  RITER  to  indicate  itera¬ 
tion  concluded 

Control  for  calling  SETUP  only  once 
Number  of  values  in  CHIT 
Number  of  values  in  CLDRT 

Number  of  points  used  to  define  pressure  vs  time  curve 

Number  of  present  trial  in  iterative  process  for  options 
1  and  2 

Number  of  values  in  UT 
Number  of  values  in  ZETAT 
Number  of  values  in  ZTA2T 
Number  of  values  in  ZETAI 
Number  of  values  in  ZTA1T 

Output  array  for  blast  characteristics  returned  by 
HYPRA 

Critical  pressure  for  Youngdahl's  method 

Critical  pressure  for  Youngdahl's  method,  metric 

Value  of  P  /P 
o  o, 
n  1 
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PH(3) 

PHP 

PI 

PO 

POI 

PP(20) 

*PPM(20) 

PPR 

PPRM 

PRFAC(102) 

*PY 

PZ(3) 

PZERO 

R 

RHO 

RHOO 

RHOT(2) 

*RM 

*ROT 

RTRIAL(3) 

SCHI 
SF  (5) 
SCHIT(31) 

*SM 


Critical  pressures  for  three  damage  levels 
Critical  pressure 
n 

Sea-level  pressure 

Ratio  of  blast  pressure  to  blast  impulse 

Value  of  pressure  at  time  TPP 

Value  of  pressure  at  time  TPP,  metric 

Peak  pressure  of  exponential  forcing  function 

Peak  pressure  of  exponential  forcing  function,  metric 

Tabulated  values  of  fitting  function  P  / P  ,  plotted 
in  Figure  4.6  n  1 

Ratio  of  critical  pressure  to  be  used  in  modified  Young- 
dahl  procedure  to  the  pressure  asymptote  for  light  damage 

Pressure  asymptotes  for  three  damage  levels 

Pressure  asymptote 

Range 

Material  density 
Sea-level  air  density 
Table  of  material  densities 
Range,  metric 

Cylinder  radius  to  thickness  ratio 

Trial  ranges  in  iteration  for  option  1;  inverse  of 
trial  pressure  for  option  2 

Value  of  s(x)  used  in  Eq.  (4.28) 

Scaling  factors  for  blast  calculations 

Tabulated  values  of  fitting  function  s(x)  plotted  in 
Figure  4.3 

Longest  distance  from  any  point  in  face  of  structure 
containing  beam  or  plate  to  farthest  free  edge 
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SMALT 


Plate  or  beam  thickness 


SN2 

SN2A 

SN2B 

SRAE 

*SSM 

STFAC (39) 

*TO 

*TM 

*TPP (20) 

U 

UT(13) 

W 

W1 

WL 

*WM 

WST 

WW 

WY 

*XX 

ZEROI 

ZETA 


Stress  parameter,  n?s(x),  in  Eq.  (4.28) 

Value  of  stress  parameter,  rrs(x),  just  after  snap- 
through 

Value  of  stress  parameter,  n?s(x),  just  before  snap- 
through 

Square  root  of  product  of  density,  area  squared,  and 
modulus  of  elasticity 

Distance  from  center  of  beam  or  plate  to  closest  free 
edge 

Tabulated  values  of  snap-through  factor  plotted  in 
Figure  4.9 

Duration  of  exponential  forcing  function 
Plate  or  beam  thickness,  metric 
Time  at  which  pressure  is  given 
P,  defined  in  Eq.  (4.9) 

Values  of  p  used  in  tabulating  data  for  snap-through 
factor  plotted  in  Figure  4.9 

Value  of  loading  parameter,  u>,  defined  in  Eq.  (4.16) 

Critical  value  of  loading  parameter,  w,  defined  in 
Eq.  (4.16),  for  light  damage 

Width  over  which  pressure  acts 

Weapon  yield,  metric 

Value  of  loading  parameter,  w,  defined  in  Eq.  (4.16),  at 
snap-through 

Critical  value  of  loading  parameter,  w,  defined  in  Eq. 
(4.16) 

Weapon  yield 
Input  dummy: 

for  beams,  width  over  which  loading  acts,  metric; 
for  plates,  aspect  ratio 

Impulse  asymptote 

Value  of  tensile  axial  load  parameter,  4,  defined  in 
Eq.  (4.19) 
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ZETAI (45) 

Values  of  C 
Figure  4.7 

used 

ZETAT (31 ) 

Values  of  C 
4.4 

used 

Z I  (3) 

Impulse  asymptote 

ZTA1T (62) 

Values  of  C 
4 . 5  and  4 . 8 

used 

ZTA2T (51) 

Values  of  C 

used 

Figure  4.6 


for  I  / I  in 

o  o, 
n  1 

for  6(c)  in  Figure 

levels 

for  6(C)  in  Figures 

for  P  /P  in 

°  °1 
n  1 


5.4  PIVUL  Input 
Card  1  (112) 

IVOPT,  vulnerability  option  code 

1- given  weapon  yield,  determine  critical  range 

2- given  duration  of  exponential  forcing  function,  deter¬ 
mine  critical  pressure  level 

3- given  weapon  yield  and  range,  determine  damage  level 

4- given  curve  of  pressure  vs  time,  determine  damage 
level 

Note:  Only  options  1  and  3  can  be  used  for  cylinders. 

Card  2  (Card  2  may  be  a  single  card  or  a  group  of  cards,  depend¬ 
ing  upon  IVOPT) 

IV0PI=1  (F12.1) 

WM,  weapon  yield*  (TJ) 

IV0PT=2  (F12.1) 

TO,  duration  of  exponential  forcing  function  (sec) 

IV0PT=3  (2F12.1) 

WM,  weapon  yield*  (TJ) 

RM,  range  (meters) 


The  program  does  not  include  ground  reflection  effects,  so  the  user 
must  adjust  the  yield  if  appropriate  to  account  for  such  effects. 
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IV0PT=4  (112,  (6F12.1)) 


NPP,  number  of  points  used  to  define  pressure  vs  time 
curve  (<20) 

TPP(N),  time  at  which  pressure  is  given  (sec) 

PPM(N),  value  of  pressure  at  time  TPP(N)  (kPa) 

The  first  card  contains  NPP  only. 

Subsequent  cards  contains  pairs  of  TPP(N)  and  PP(N), 
three  pairs  to  a  card.  A  total  of  NPP  pairs  must  be 
provided.  The  first  time  must  be  0.0.  The  last  pressure 
must  be  0.0. 

Card  3  (F12.1) 

PY,  ratio  of  critical  pressure  to  be  used  in  modified  Young- 
dahl  procedure  to  P  ,  the  pressure  asymptote  for  light 
°1 

damage.  It  is  recommended  that  PY  be  assigned  a  value 
of  0.75. 

If  PY>1.0,  the  modified  Youngdahl  procedure  is 

not  used.  PY  must  be  positive.  For  cylinders,  PY  is 

not  used,  and  any  value  may  be  put  on  the  input  card. 

Card  4  (3112) 

IMAT,  material  code 

1- 2024T4  Aluminum 

2- 7075T6  Aluminum 

IBPC ,  structure  code 

1- beam 

2- panel 

3- cylinder 

ISL,  damage  level  code  (required  only  for  IV0PT=1  or  2) 

1- light  damage 

2- moderate  damage 

3- heavy  damage 
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The  following  card(s)  define  the  structure  and  some  aspects  of  the 
loading.  They  depend  upon  the  code  IBPC,  and  will  be  given  for  each  of 
the  three  possible  values  of  IBPC. 

IBPC=1 

Card  5  (5F12.1) 

BLM,  beam  length  (meters) 

CLDR,  length  to  radius  of  curvature  ratio  (0.0<CLDR_^0.8) 

XX,  width  over  which  loading  acts  (meters) 

SSM,  distance  from  center  of  beam  to  closest  free  edge  (meters) 
(omit  for  IV0PT=4) 

SM,  longest  distance  from  any  point  in  face  of  structure 
containing  beam  to  farthest  free  edge  (meters) 

(omit  for  IV0PT=4) 

Card  6  (4F12.1) 

ATOTM,  cross-section  area  (sq  cm) 

4 

CNTIM,  cross-section  area  moment  of  inertia  (cm  ) 

CNTRM,  distance  from  neutral  axis  to  outer  fiber  on  loaded 
side  (cm) 

TM,  beam  thickness  (cm) 

This  completes  the  input  for  IBPC=1. 

IBPC=2 

Card  5  (5F12.1) 

BLM,  length  of  shorter  side  of  panel  (meters) 

CLDR,  ratio  of  BLM  to  radius  of  curvature  in  direction  corre¬ 
sponding  to  BLM  (0.0<CLDRf0. 8) 

AR,  aspect  ratio  of  plate,  ratio  of  length  of  longer  side 
to  length  of  shorter  side 

SSM,  distance  from  center  of  plate  to  closest  free  edge 
(meters)  (omit  for  IV0PT=4) 

SM,  longest  distance  from  any  point  in  face  of  structure  con¬ 
taining  plate  to  farthest  free  edge  (meters) (omit  for 
IV0PT=4) 


Card  6  (F12.1) 


TM,  plate  thickness  (cm) 


This  completes  the  input  for  IBPC=2. 


IBPC=3 


Card  5  (2F12.1) 

ROT,  radius  to  thickness  ratio  (5(KROT<)00) 
CLDD,  length  to  diameter  ratio  (1<CLDD<_3) 

This  completes  the  input  for  IBPC=3. 


5 . 5  Program  Operation  and  Output 

The  NOVA  program  is  written  in  FORTRAN-IV  and  has  been  run  on  the 
Control  Data  Corporation  6600  scientific  computer  under  the  SCOPE 
(Supervisory  Control  of  Program  Execution)  operating  system,  version 
3.4.3.  About  70,000  cells  of  central  memory  are  required  to  compile 
PIVUL  and  40,000  cells  of  central  memory  are  required  for  program  load¬ 
ing  and  execution.  Computation  time  is  on  the  order  of  0.2  seconds. 

Although  the  program  output  is  largely  self-explanatory,  the  normal 
output  is  described  in  detail  in  Table  5.1.  The  corresponding  program 
variable  is  given  parenthetically. 

In  addition,  certain  error  situations  can  develop  which  will 
generate  error  messages,  which  are  described  in  Table  5.2.  In  this 
connection,  it  should  be  noted  that  the  user  is  responsible  for  observ¬ 
ing  the  limits  on  input  variables  indicated  in  Subsection  5.4.  Gener¬ 
ally,  erroneous  input  will  not  generate  an  error  message.  The  input,  or 
information  directly  devised  from  input,  is  printed  out,  however,  so 
that  the  user  should  easily  be  able  to  identify  that  his  problem  is  due 
to  incorrect  input. 
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TABLE  5.1 
NORMAL  OUTPUT 


VULNERABILITY  OPTION  i  (IVOPT) 


For 


For 


For 


For 


IV0PT=1 , 

YIELD  =  X.XXXXXX  TJ  (WM) 

IVOPT=2, 

DURATION  =  X.XXXXXX  SEC  (TO) 

IVOPT=3, 

YIELD  =  X.XXXXXX  TJ  AT  RANGE  OF  X.XXXXXX  METERS  (WM,  RM) 
IV0PT=4, 

TIME  PRESSURE 

(SEC)  (KPA) 

X.XXXXXX  X.XXXXXX 


i 

i 

i 


i 

I 

I 


X.XXXXXX  X.XXXXXX 


YOUNGDAHLS  METHOD,  PCR/PO  =  X.XXXXXX  (PY,  output  only  if  PY<1.0) 
LIGHT  DAMAGE  or  MODERATE  DAMAGE  or  HEAVY  DAMAGE 
2024-T4  ALUMINUM  or  7075-T6  ALUMINUM 
BEAM  or  PANEL  or  CYLINDER 


For  beams  (IBPC=1) , 

LENGTH  =  X.XXXXXX  METERS  (BLM) 

L/R  =  X.XXXXXX  (CLDR) 

LOADING  WIDTH  =  X.XXXXXX  METERS  (XX) 

*  SHORTEST  DISTANCE  TO  FREE  EDGE  =  X.XXXXXX  METERS  (SSM) 

*  LONGEST  DISTANCE  FROM  ANY  POINT  TO  FREE  EDGE  =  X.XXXXXX  METERS  (SM) 
CROSS-SECTION  DATA 

AREA  =  X.XXXXXX  SQ  CM  (ATOTM) 

MOMENT  OF  INERTIA  =  X.XXXXXX  CM**4  (CNTIM) 

DISTANCE  TO  OUTSIDE  FIBER  ON  LOADED  SIDE  =  X.XXXXXX  CM  (CNTRM) 
THICKNESS  =  X.XXXXXX  CM(TM) 


*Omitted  for  IVOPT= 


4. 
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TABLE  5.1  (Concluded) 


For  panels  (IBPC=2) 

LENGTH  =  X.XXXXXX  METERS  (BLM) 

L/R  =  X.XXXXXX  (CLDR) 

ASPECT  RATIO  =  X.XXXXXX  (AR) 

*  SHORTEST  DISTANCE  TO  FREE  EDGE  =  X.XXXXXX  METERS  (SSM) 

*  LONGEST  DISTANCE  FROM  ANY  POINT  TO  FREE  EDGE  =  X.XXXXXX  METERS  (SM) 
THICKNESS  *  X.XXXXXX  CM  (TM) 

For  cylinders  (IBPC=3) 

R/T  =  X.XXXXXX  (ROT) 

L/D  =  X.XXXXXX  (CLDD) 

PCR  FOR  YOUNGDAHLS  METHOD  =  X.XXXXXX  KPA  (PCRM,  output  only  if  PY<1.0) 
For  IVOPT=l, 

CRITICAL  RANGE  =  X.XXXXXX  METERS  (RM) 

FOR  IVOPT=2 , 

CRITICAL  FREE-FIELD  OVERPRESSURE  =  X.XXXXXX  KPA  (PPRM) 

For  IVOPT=3  or  4, 


ZERO 

LIGHT 

MODERATE 

HEAVY 

RESPONSE 

DAMAGE 

DAMAGE 

DAMAGE 

Y 

Y 

Y 

Y 

X 


* 

Omitted  for  IV0PT*4* 
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TABLE  5.2 
ERROR  MESSAGES 


CHI  =  X.XXXXXX  CHECK  BEAM  CROSS-SECTION  DATA 

This  message  is  printed  out  by  BEAM  if  CHIS1  is  not  possible, 
the  beam  cross-section  data  must  be  in  error.  Program  stops. 

C  VS  R  CURVE  SLOPE  REVERSAL 

This  message  is  printed  out  by  RITER  if  the  response  criterion, 
CRIT,  has  decreased  with  decreasing  range  (or  increasing  pressure). 
This  violates  the  assumption  of  RITER  that  the  criterion  increases 
monotonically  with  increasing  pressure.  Program  stops. 

ERROR  STOP  -  CRITICAL  PRESSURE  FOR  YOUNGDAHL  IS  GREATER  THAN  PEAK 
APPLIED  PRESSURE 

For  option  3  or  4,  the  peak  pressure  in  the  forcing  function  is 
smaller  than  PCR,  thus  making  Youngdahl  integrations  invalid. 

IVULT  prints  out  this  message  and  program  stops. 

ERROR  STOP  -  DAMAGE  LEVEL  iiiiii  NOT  ALLOWED 

The  damage  level  code  is  less  than  1  or  greater  than  3.  PIVUL 
prints  out  this  message  and  program  stops. 

ERROR  STOP  -  DIFFRACTION  LOADING  PERIOD  TOO  LONG  -  CHECK  SM  AND  SSM 
This  message  is  printed  out  by  SIDE  if  the  diffractive  clearing 
time  is  greater  than  the  forcing  function  duration.  User  should 
check  values  of  input  variables  SM  and  SSM.  Program  stops. 

FIRST  ARGUMENT  IS  OUTSIDE  TABLE,  X=X.XXXXXX 
XT  =  X.XXXXXX . X.XXXXXX 

This  message  is  printed  out  by  INTERP  if  the  first  independent 
variable,  X,  falls  outside  the  range  of  tabulated  data.  The 
value  of  the  first  independent  variable  and  the  corresponding 
table  values  are  printed  out.  Program  stops. 

RANGE  ITERATION  IS  NOT  CONVERGING 

To  avoid  a  possible  iteration  loop  in  searching  for  the  critical 
range  or  pressure,  the  number  of  trials  allowed  is  restricted  to 
twelve.  This  message  is  printed  out  by  RITER  to  indicate  that 
adequate  convergence  has  not  been  achieved  in  twelve  trials. 

Program  stops. 

SECOND  ARGUMENT  IS  OUTSIDE  TABLE,  Y=X.XXXXXX 
YT  =  X.XXXXXX . X.XXXXXX 

This  message  is  printed  out  by  INTERP  if  the  second  independent 
variable,  Y,  falls  outside  the  range  of  tabulated  data.  The  value 
of  the  second  independent  variable  and  the  corresponding  table 
values  are  printed  out.  Program  stops. 
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TABLE  5.2  (Concluded) 


THIRD  ARGUMENT  IS  OUTSIDE  TABLE,  Z  =  X.XXXXXX 
ZT  =  X.XXXXXX . X.XXXXXX 

This  message  is  printed  out  by  INTERP  if  the  third  independent 
variable,  Z,  falls  outside  the  range  of  tabulated  data.  The 
value  of  the  third  independent  variable  and  the  corresponding 
table  values  are  printed  out.  Program  stops. 

VULNERABILITY  OPTION  i  NOT  ALLOWED  FOR  CYLINDERS 

Vulnerability  option  2  or  4  has  been  specified  for  a  cylinder. 
PIVUL  prints  out  this  message  and  program  stops. 

X- VALUE  X.XXXXX  OUTSIDE  TABLE.  Y-VALUE  RETURNED 

=  X.XXXXXX  X-TABLE  =  X.XXXXXX . X.XXXXXX 

This  message  is  pointed  out  by  INTI  if  the  independent  variable, 

X,  falls  outside  the  range  of  tabulated  data.  The  value  of  the 
independent  variable  and  the  table  values  are  printed  out.  The 
dependent  variable,  Y,  will  be  the  minimum  or  maximum  table  value, 
depending  upon  whether  X  falls  below  or  above  the  tabulated  data. 
Program  stops. 
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5 . 6  Example  Problem 

The  structure  for  the  example  problem  is  a  curved  panel.  The 
required  dimensions  and  material  are: 

Thickness  =  0.08  cm 

Length  of  short  side  =0.2  meters 

Length  to  radius  of  curvature  ratio  =  0.05 

Aspect  ratio  =  1.3 

Material  =  2024T4  A1 

Distance  from  center  of  plate  to  closest  free  edge  =  1  meter 

Longest  distance  from  any  point  in  face  of  structure  contain¬ 
ing  plate  to  farthest  free  edge  =  3  meters 

It  is  desired  to  determine  the  level  of  response  of  this  structure  at  a 
range  of  300  meters  from  a  4.2  TJ  explosion  (the  value  4.2  is  an  effec¬ 
tive  yield  which  includes  enhancement  by  ground  effects).  The  modified 
Youngdahl  approach  is  to  be  used,  with  a  ratio  of  critical  pressure 
to  P  of  0.75. 

°1 

The  input  data  deck  is  as  follows: 

3 

4.2  300. 

0.75 

1  2 

0.2  0.05  1.3  1.0  3.0 

0.08 
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The  resulting  output  is  reproduced  below. 


Vllt-Ntn  Ab  1 L  1  I  Y  OHIiO*  S 

tltLli  S  . A2U 0 00 t ♦ U 1  rj  A(  KAubt  lit-  .iUOOUUt  +  OJ  l,'tltK3 

YOUNGiUAhLS  Mtlhdd,  kCk/HO  =  ,7bu0uut  +  u0 

«;02M-IO  ALUMINUM 

P  A  .Ni  E  U 

s  .iOOOoot+VU  MtltKS  u / K  s  .bOUOUOc-ul 

ASPfcCl  KA | 1U  =  . 1 4dUUUt+Ul 

tUbt  =  .10oUUUfc.+Ul  Mtltwb 
Hulm?  1  U  H<  t  fc  t  Dbt  =  .i0vuGufc*01  MtltKb 
CM 

s  .iiS4bbt  +  Uei  AKA 

Li1**1  wUOfcKATt  HLAVT 

OAMAGl  U  A  M  A  C"  t  u  A  r,  A  G  ' 

'  t  r 

A 


SMuHTtbT  OlSTANCt  1U  t-Kfct 
LUNbtbl  DialANCt  f-'Huw  AmY 
TnlCKAtbb  =  .ojuuout-01 

HCH  Fuh  Y  UUiMbDAntb  htltlUU 

JfcKu 

HtbKUn8t 

Y 


The  output  shows  that  the  structure  experiences  an  amount  of  damage 
between  moderate  and  heavy  for  the  specified  conditions. 
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SECTION  6 

CONCLUSIONS  AND  RECOMMENDATIONS 


This  report  has  presented  the  methodology  which  led  to  the  develop¬ 
ment  of  a  vulnerability  code  for  utilization  in  rapid  damage  assessment 
of  simple  structural  elements  subjected  to  a  nuclear  blast  environment. 
In  the  development  of  the  methodology,  structural  models  typifying  a 
variety  of  representative  straight  and  curved  beam,  flat  and  curved 
plate,  and  cylinder  configurations  were  analytically  loaded  with  simu¬ 
lated  nuclear  blast  overpressure  time-histories  and  light,  moderate,  and 
heavy  damage  levels  of  response  were  defined  by  using  pressure-impulse 
analyses.  The  structural  model  support  constraints  considered  were 
fixed-fixed  end  conditions  for  the  beams,  plates,  and  cylinders;  the 
materials  of  construction  considered  were  both  2024-T4  and  7075-T6 
aluminum  alloys.  During  the  course  of  the  investigation,  several 
general  conclusions  and  recommendations  were  arrived  at  and  are  pre¬ 
sented  in  the  paragraphs  which  follow. 

The  moment  gradient  which  exists  at  the  ends  of  a  thin  fixed-fixed 
beam  is  extremely  steep,  due  to  the  assumed  boundary  condition  of  a 
perfectly  rigid  support.  This  steep  moment  gradient  in  turn  results  in 
localized,  extremely  severe  strains  at  the  ends  of  the  beam.  Response 
criteria  based  on  these  highly  localized  strains  may  be  too  severe  in 
terms  of  the  improbability  of  attaining  a  perfectly  rigid  support  in 
real  construction.  Edge  fastenings,  which  are  not  modeled,  may  be  of 
importance. 

The  use  of  static  elastic  analysis  solutions  as  a  guide  in  selec¬ 
ting  appropriate  fitting  functions  for  the  P-I  data  generated  proved  to 
be  very  helpful.  The  particular  success  achieved  in  applying  the  method 
to  the  large  amount  of  2024  straight  beam  data  generated  resulted  in  a 
fairly  high  degree  of  confidence  in  the  resulting  2024  straight  beam 
fitting  functions.  Somewhat  less  confidence  is  assigned  to  the  2024 
curved  beam  fitting  functions  due  to  the  added  complexity  posed  by  the 
curved  beam  snap-through  problem. 

Since  a  lesser  amount  of  P-I  data  were  generated  for  7075  struc¬ 
tures  than  for  2024  structures,  it  was  necessary  to  employ  conversion 
factors  to  relate  the  selected  7075  fitting  functions  to  the  previously 
selected  2024  fitting  functions.  Consequently,  the  7075  fitting  func¬ 
tions  are  considered  to  be  not  quite  as  good  as  the  2024  fitting  func¬ 
tions.  However,  if  more  7075  P-I  data  were  available,  static  elastic 
analysis  could  again  be  utilized  as  a  guide  in  selecting  7075  fitting 
functions  which  would  result  in  a  higher  degree  of  confidence. 
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The  vulnerability  assessment  code  which  has  been  developed  herein 
is  satisfactory  for  providing  rapid,  approximate  answers  as  to  the 
relative  hardness  of  simple  structural  elements  subjected  to  nuclear 
blast  overpressure  loadings.  Furthermore,  the  code  should  prove  to  be  a 
valuable  tool  when  used  either  in  the  preliminary  designs  of  future 
military  weapon  systems  or  in  the  analysis  of  presently  existing  systems. 

Typical  of  any  analytical  tool  development,  the  vulnerability  code 
presented  herein  requires  experimental  verification.  It  is  recommended 
that  experimental  verification  begin  with  the  2024  straight  beam.  If 
such  experimental  verification  is  successful,  experiments  should  be  run 
for  correlation  with  the  2024  curved  beam  and  2024  plate  analyses. 

If  the  experimental  verification  is  successful,  further  analytical 
efforts  should  be  undertaken  to  obtain  more  accurate  fitting  functions 
for  the  7075  aluminum  elements  and  to  develop  data  for  other  materials. 
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APPENDIX  A 

VARIATION  IN  P-I  CURVE  WITH  FORCING  FUNCTION 


A  pressure-impulse  curve  is  constructed  for  a  given  forcing  func¬ 
tion  time-history,  and  a  P-I  curve  should  be  used  cautiously  for  other 
forcing  functions.  This  appendix  addresses  the  question  of  how  much  the 
P-I  curve  varies  with  forcing  function  for  a  simple  class  of  forcing 
functions  typical  of  blast-type  loadings. 

The  forcing  function  used  for  this  purpose  is  shown  in  Figure  A.l, 
and  consists  simply  of  the  superposition  of  two  triangular  forcing  func¬ 
tions  of  different  durations.  The  amplitude  of  the  longer  duration 
triangle  is  one-half  the  amplitude  of  the  total  forcing  function  at  time 
zero.  The  parameter  which  is  used  to  obtain  different  forcing  functions 

is  the  duration  ratio,  t,/t  . 

1  o 

The  structure  for  which  the  forcing  functions  were  applied  was  a 
2024-T4  aluminum  I-beam,  with  a  length  to  height  ratio  of  35  and  area 
ratios  A^ /A^  =  0.9  and  A^/A^  =  1.0.  The  light  damage  case  was  examined. 

The  resulting  P-I  curves  are  given  in  Figure  A. 2  for  duration 
ratios  of  0.01,  0.5,  and  1.0.  The  pressure  and  impulse  are  normalized 
to  the  pressure  and  impulse  asymptotes,  respectively.  The  range  of 
duration  ratios  covers  the  extremes:  The  value  of  0.01  corresponds  to  a 
very  sharp  spike  at  the  beginning,  while  the  value  of  1.0  results  in 
just  a  simple  triangular  forcing  function. 

The  variation  between  the  two  extreme  curves  is  seen  to  be  signi¬ 
ficant.  To  quantify  the  significance,  suppose  an  analyst  has  a  pressure¬ 
time  curve  of  given  duration  and  wants  to  know  the  magnitude  of  the 
pressure  required  to  produce  damage.  A  given  duration  defines  a  straight 
line  through  the  origin  of  the  P-I  plane.  Hence,  the  ratio  of  the 
pressure  from  the  two  extreme  curves  along  a  line  through  the  origin  is 
the  desired  measure  of  the  significance  of  the  variation  in  the  P-I 
curves.  The  largest  value  of  this  ratio  occurs  along  the  line  defined 
by  the  normalized  pressure  and  the  normalized  impulse  being  equal  (a  45° 
line  if  the  horizontal  and  vertical  scales  are  the  same,  as  in  Figure 
A. 2),  and  is  approximately  1.4.  Hence,  if  the  analyst  had  a  P-I  curve 
based  on  a  duration  ratio  of  0.01  available  and  was  interested  in  a 
simple  triangular  forcing  function  (duration  ratio  =  1),  he  would  over¬ 
predict  the  pressure  required  by  40%  at  most.  If  he  had  a  P-I  curve 
midway  between  the  extreme  curves  on  Figure  A. 3  he  could  not  be  in  error 
by  more  than  20%. 

It  should  be  remembered  that  the  above  errors  result  when  consi¬ 
dering  only  forcing  functions  of  the  class  depicted  by  Figure  A.l.  For 
other  forcing  functions,  such  as  a  ramp  or  a  half  sine,  the  error  could 
be  significantly  greater. 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet  and  return  it  to  Director,  US  Army  Ballistic  Research 
Laboratory,  ARRADCOM,  ATTN:  DRDAR-TSB,  Aberdeen  Proving  Ground, 
Maryland  21005.  Your  comments  will  provide  us  with  information 
for  improving  future  reports. 

1.  BRL  Report  Number^ _ _ _ _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) _ _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 


Name: 

Telephone  Number: 
Organization  Address: 


